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The  dynamics  of  the  large  scale  structures/vorticlty  concentrations  ob¬ 


J 


served  in  two-dimensional  transitional  and  turbulent  mixing  layers  are 
described  in  terms  of  a  nonlinear  invlscid  twcjjdimensional  instability 
model.  According  to  the  model; the  formation  and  repetitive  pairing  of 
the  structures  reflect  a  deterministic  cyclic  mechanism  wherein  the  non¬ 
linear  evolution  of  the  flnitelamplitude  instability,which  dominates  the 
flow  at  any  one  stage,  intrinsically  provides  first  excitation  of,  and 
then  energy  transfer  to,  the  bound  1  /id  subharmonic.  The  key  to  identify-  j 

ing  and  modeling  the  mechanism  resides  in  a  realistic  description  of  the  con¬ 
current  changes  in  the  mean  flow  and  in  the  structure  of  the  wave  motion 
during  the  amplification,  saturation  and  subsequent  modulation  of  the  indi¬ 
vidual  nonlinear  instability.  Models  ba^ed  on  weakly  nonlinear  stability 
theory  prove  inadequate  to  this  task.  The  cause  of  the  failure  is  Identi¬ 
fied,  and  attendant  criteria  for  estimating  the  applicability  of  the  theory 
to  the  present  as  well  as  other  problems  are  indicated.  By  contrast,  several 
salient  features  of  transitional  and  turbulent  mixing  layer  development  are 
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reproduced  by  the  repeated  application  of  a  nonlinear  model, which  describes  the 
life  cycle  of  the  Individual  structure  only  in  terms  of  the  interactive  evolu¬ 
tion  of  the  mean  flow  and  of  the  wave  motions  associated  with  an  appropriate 


finite  amplitude  fundamental  together  with  its  bbund  second  harmonic 


Some 


views  on  possible  external  controls  influencing  the  development  of  the  struc 


tures  and  of  the  attendant  mean  flow  are  Inferred  on  that  basis. 


1.  INTRODUCTION 

Recent  experiments  by  Wygnanski  et  al  (1979a)  and  by  Browand  and  Troutt  (1980) 
clearly  demonstrate  that  the  quasi  two  dimensional  large  scale  coherent  struc¬ 
tures/  vorti  city  concentrations  first  observed  by  Brown  and  Roshko  (1974)  per¬ 
sist  In  turbulent  mixing  layers.  The  structures,  with  principal  alignment 
In  the  spanwise  direction,  rapidly  approach  an  asymptotic  state  where  their 
mean  transversal  dimensions  as  well  as  spanwise  correlation  length  become 
functions  only  of  the  local  mixing  layer  thickness  6^.  Their  mean  distance  of 
travel  between  pairings  then  depends  only  on  the  ratio  (6^/X)  =  6w(U-j  +  Up)/ 

(U-j  -  U2)— where  Uj  and  U2  denote  the  velocities  of  the  two  streams  partaking 
in  the  mixing* -so  that,  for  a  fixed  velocity  ratio  X,  the  length  of  the  pairing 
process  also  scales  with  the  local  thickness  6^.  Thus,  the  repetitive  occur¬ 
rence  of  a  deterministic  self-preserving  mechanism  of  intrinsic  subharmonic 

\ 

generation  is  strongly  suggested.  The  view  is  supported  by  the  observation 
that  initially  forced  mixing  layers  (Wygnanski  et  al  1979b,  Ho  and  Huang  1980) 
exhibit  a  rapid  return  to  the  behavior  of  unforced  layers  following  the  initial 
formation  of  a  single  large  structure  by  the  simultaneous  agglomeration  of 
several  smaller  ones. 

Whereas  the  pairing  of  structures/vortices  is  largely  responsible  for  the  turbu¬ 
lent  mixing  layer  growth  (Winant  and  Browand  1974,  Roshko  1976),  a 
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quantitative  understanding  of  the  attendant  dynamics  Is  of  distinct  Interest 
from  the  dual  viewpoint  of  developing  flow  models  based  on  first  principles  and 
of  establishing  a  rational  basis  for  possible  flow  control  concepts.  Available 
models  of  the  pairing  process  rely  on  numerical  analyses  of  the  unsteady 
twodlmenslonal  flows  which  evolve  either  from  an  Initial  parallel  "base  shear 
flow"  perturbed  by  one  or  more  excited  elgenmodes  of  arbitrarily  prescribed  small 
amplitude  end  relative  phase  (e.g.  Zabusky)and  Deem  1971,  Patnalk  et  al  1976), 
or  from  an  assembly  of  discrete  line  vortices  Initially  arranged  In  cyclic 
arrays  with  arbitrarily  prescribed  small  lateral  displacements  and/or  offsets 
(e.g.,  Acton  1976),  Aref  and  Slggla  1980)  .  In  either  case  solutions  are 
obtained  subject  to  boundary  conditions  which  assume  periodicity  of  all  flow 
variables  in  the  streamwise  direction  over  the  length  of  the  computational 
domain.  Many  detailed  features  experimentally  observed  during  the  nonlinear 
phases  of  transition  In  selected  free  shear  flows  have  successfully  been  repro¬ 
duced  by  these  models  in  conjunction  with  Initial  conditions  which  realistically 
reflect  those  encountered  in  the  experiments.  However,  the  comparisons  have 
typically  been  restricted  to  limited  Intervals  of  space  in  the  experiments  and 
corresponding  time  in  the  theoretical  predictions.  With  few  exceptions,  e.g., 
the  calculation  reported  by  Aref  and  Siggia  (1980),  the  onset  of  a  cyclic  self¬ 
preserving  process  of  1/2  subharmor.lc  generation  and  amplification,  such  as 
Indicated  by  the  observations  of  coherent  structures  in  mixing  layers,  has  not 
been  demonstrated  analytically.  Since  the  numerical  solutions  are  deterministic 
and  uniquely  defined  by  the  assumed  initial  conditions  and  the  non  controversial 
boundary  conditions,  it  is  difficult  to  understand  how  they  could  predict  self¬ 
preserving  vortex  pairing  over  distances  which  scale  with  the  local  layer 
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thickness  In  unforced  mixing  layers,  and  rapidly  recover  the  same  self- 
preserving  process  In  forced  mixing  layers,  which  evolve  from  obviously  differ¬ 
ent  Initial  conditions.  Thus,  the  Issue  of  understanding  the  dynamics  of 
the  scale  change  undergone  by  the  structures  as  they  convect  downstream  Is 
not  fully  settled. 

The  observed  rapid  recovery  of  self -preserving  vortex  pairing  over  lengths  which 
scale  with  the  local  mixing  layer  thickness,  Independently  of  the  Initial  and 
external  disturbance  spectrum,  suggests  that  the  operative  mechanism  may  be 
one  where  I)  the  excitation  of  the  1/2  subharmonic  Is  intrinsic  to  the  non¬ 
linear  evolution  of  the  finite  amplitude  Instabilities  which  control  the  forma¬ 
tion  of  the  individual  structures,  and  2)  the  subsequent  amplification  of  the 
subharmonic  drains  energy  from  the  preexisting  instability  so  that  memory  of 
the  specific  Initial  conditions  Is  largely  lost  after  one  or  a  few  structure 
life  cycles.  This  is  the  proposition  elaborated  in  this  paper.  Its  demon¬ 
stration  clearly  depends  on  a  realistic  description  of  the  history  and  struc¬ 
ture  of  nonlinear  wave  motions  associated  with  Individual  finite  amplitude 
instabilities  and  on  a  study  of  self-induced  modulations  which  result  In 
interacting  pairs  of  such  instabilities  having  frequency  ratios  of  about  2. 

The  paper  devotes  much  attention  to  the  analysis  of  these  specific  problems. 

In  order  to  avoid  the  large  labor  and  cost  entailed  by  simulations  of  spa¬ 
tially  developing  flows,  where  the  streamwise  extent  of  the  computational 
domain  must  encompass  several  vortex  pairing  lengths,  our  investigations  address 
the  problem  of  temporally  developing  mixing  layers,  where  a  computational 
domain  with  streamwise  extent  restricted  to  two  eddy  sizes  allows  analysis  of 
the  vortex  pairing  process  and  assessment  of  its  self  preserving  nature.  In 
doing  so  we  recognize  that  the  time  evolution  of  the  unsteady  shear  layer. 
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statistically  homogeneous  In  the  streamwlse  direction,  and  the  space  evolution 
of  the  laboratory  mixing  layer,  statistically  stationary  In  time,  are  not 
strictly  equivalent, si  nee  one  flow  cannot  be  reduced  to  the  other  by  a 
Galilean  transformation.  Accordingly,  we  employ  due  caution  In  transposing 
results  from  one  situation  to  the  other,  especially  for  what  concerns  the 
Interactions  among  Instability  pairs,  which  are  known  to  be  dispersive  In  the 
spatial  amplification  case  (Mlchalke  1965)  but  non-dlsperslve  In  the  temporal 
amplification  case  (Mlchalke  1964).  The  relevance  of  this  caveat  to  the  study 
of  vortex  pairing  Is  evidenced  by  Petersen's  (1978)  experimental  finding 
that"...  the  first  pairing  occurred  where  the  mixing  layer  became  sufficiently 
thick,  owing  to  nonlinearities  and  entrainment,  that  a  subharmonic  wave  had 
a  phase  velocity  that  matched  the  convection  speed  of  the  vortex  ring." 

Our  inquiries  Into  the  structure  of  the  wave  motion  associated  with  instabilities 
of  small  but  finite  amplitude,  growing  at  a  finite  rate,  largely  address  the 
central  tenet  of  weakly  nonlinear  theories,  viz.,  that  the  eigenfunctions 
given  by  linear  theory  continue  to  yield  an  adequate  representation  of  the 
motion  in  the  nonlinear  case.  Since  experiments  (e.g.,  Mlksad  1972)  clearly 
show  that  significant  variations  of  mean  flow  and  maximum  shear  accompany  the 
growth  of  instabilities  even  at  modest  amplitudes,  we  relinquish  the  classical 
assumption  of  invariant  eigenfunction  shape  (Stuart  1960)  in  favor  of  the 
Meksyn  ard  Stuart  approach  (1951),  which  we  rederive  in  the  context  of  multiple 
scales  methods  (e.g.,  Nayfeh  1973)  together  with  a  self-consistent  equation 
for  the  time  rate  of  change  of  the  disturbance  energy.  We  then  assess  the  pre¬ 
dictions  of  the  extended  Meksyn  and  Stuart  model  by  comparing  them  with  Miksad's 
(1972)  hot  wire  measurements  in  the  nonlinear  stages  of  free  shear  layer  transi¬ 
tion  as  well  as  with  the  predictions  of  a  nonlinear  model  largely  patterned 
after  those  of  Zabusky  and  Deem  (1971)  and  Patnaik  et  al  (1976).  Unfortunately, 
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the  assessment  proves  negative,  so  that  all  our  subsequent  considerations 
necessarily  rely  on  the  nonlinear  model.  On  the  positive  side,  the  reason 
for  the  failure  of  the  Meksyn  and  Stuart  approach  Is  Identified,  and 
attendant  criteria  for  estimating  the  applicability  of  that  model  In  other 
nonlinear  Instability  problems  are  Indicated. 

Our  Inquiries  Into  the  mechanism  of  subharmonic  generation  and  amplification 
largely  stem  from  the  predictions  of  the  nonlinear  model  for  the  concurrent 
evolutions  of  the  mean  flow  and  of  the  wave  motion!  associated  with  a  finite 
amplitude  instability  undergoing  saturation  and  subsequent  self-induced  modu¬ 
lation  of  both  amplitude  and  structure.  Whereas  the  frequency  of  the  modula¬ 
tion  is  close  to  that  of  the  subharmonic,  and  the  phase  velocity  of  the  sub¬ 
harmonic  compatible  with  the  prevailing  mean  flow  structure  matches  the  con¬ 
vection  speed  of  the  vortlcity  concentrations  generated  by  the  saturated 
instability,  we  find  here  conditions  intrinsically  conducive  to  the  initiation 
of  vortex  pairing.  A  calculation  employing  the  nonlinear  model  to  describe 
the  development  of  the  unsteady  shear  layer  which  evolves  from  such  Initial 
conditions  (including  the  attendant  Initial  structures  of  the  fundamental  and 
subharmonic  wave  motions)  allows  us  then  to  follow  the  pairing  process  to 
subharmonic  saturation,  and  to  demonstrate  the  cyclic  nature  of  that  process 
by  the  favorable  comparison  of  the  mean  flow  and  wave  motion  prevailing  at 
the  beginning  and  at  the  end  of  the  cycle. 

In  the  main  body  of  the  paper,  the  elaboration  of  our  inquiries,  and  the  pre¬ 
sentation  of  the  attendant  results,  are  prefaced  by  a  derivation  of  the  relevant 
analytical  models  in  the  context  of  selected  experimental  data,  which  provide 
guidelines  for  that  derivation  (Section  2),  as  well  as  benchmarks  for  the 
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subsequent  assessments  (Sections  4  and  5).  Numerical  procedures  used  In  the 
implementation  of  the  models  are  presented  in  Section  3. 

2.  FORMULATION  OF  THE  MODELS 

Based  on  the  experimental  observation  (Roshko  1976)  that  the  large  organized 
structures  of  mixing  layers  are  not  affected  by  the  small  scale  turbulence 
appearing  at  the  higher  Reynolds  numbers. we  take  the  threefold  view  that: 

1)  essentially  the  same  cyclic  nonlinear  mechanism  governs  the  formation  and 
the  pairing  of  structures  throughout  the  nonlinear  transitional  and  turbulent 
regimes;  2)  In  each  cycle  the  mechanism  is  only  associated  with  a  few  discrete 
large  scale  dominant  Invlscid  twodimensional  Instabilities, whose  energy  first 
grows  and  then  decays  along  the  lines  Indicated  by  measurements  in  transi¬ 
tional  flows;  3)  the  high  Reynolds  number  small  scale  turbulence,  bound  to  ' 

the  large  scale  structures  and  generated  by  a  hierarchy  of  secondary  insta¬ 
bilities  intrinsic  to  the  flow  patterns  associated  therewith,  plays  a  negli¬ 
gible  role  in  the  dynamics  of  the  individual  cycle,  and  merely  provides  the 
conduit  for  cascading  the  residual  energy  of  large  scale  instabilities  dominant 
in  previous  cycles.  In  support  of  these  views  we  submit  that  the  predominance 
of  a  two-dimensional  mechanism  is  not  precluded  by  the  experimental  finding 
(Browand  and  Troutt  1980)  that  the  structures  are  not  strictly  spanwlse  and 
two-dimensional  but,  instead,  may  be  skewed  and  branched.  The  observed  ratio 
between  the  spanwise  correlation  length  of  the  structures  and  the  local 
mixing  layer  thickness  approaches  a  lower  bound  of  about  3  for  20%  corre¬ 
lation  coefficient  in  flows  characterized  by  a  velocity  ratio  A  »  1  and  by  a  suffi 
ciently  high  Reynolds  number  (Uj-U^/v.  In  these  flows  the  streamwise  separa¬ 
tion  between  structures,  as  well  as  their  travel  distance  between  pairings. 
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are  of  the  order  of  3  to  46u;  thus,  the  wavelength  characteristic  of  the  ob¬ 
served  spanwise  skewness  Is  distinctly  larger  than  the  lengths  characteristic 
of  the  streamwise  instabilities  controlling  the  evolution  of  the  structures.* 
Under  those  conditions  only  a  linear,  zero  average,  spanwise  modulation  of 
the  |'  instabilities  and  their  nonlinear  dynamics  should  accompany, and  manifest 
itself  1n»the  skewness;  accordingly,  consideration  of  the  nonlinear  processes 
in  a  spanwise-averaged  twodimensional  approximation  is  indicated.  In  that 
connection  we  also  remark  that  a  secondary  Taylor-type  instability— |intrinsic 
to,  and  bound!  with, the  flow  pattern  relative  to  the  evolving  vorticity  concen- 
trations/vorticity  braids  characteristic  of  the  structures-^is  presumably 
responsible  for  the  observed  skewness.  Pending  a  quantitative  analysis  of 
this  instability  and  its  interactions  /with  the  streamwise  instabilities,  we 
note  that  the  experimentally  observed  invariant  scaling  of  the  spanwise  correla¬ 
tion  length  with  the  layer  thickness  <5^  indicates  synchronous  cycles  of  period 
doubling  and  attendant  energy  transfer  in  wavenumber  space  for  the  dominant 
streamwise  and  spanwise  instabilities.  Since  period  doubling  implies  at  least 
a  half  revolution  of  fluid  elements  within  the  developing  subharmonic  cat’s  eye, 
and  since  spanwise  and  streamwise  r.m.s..  velocities  are  at  best  comparable,  the 
shorter  wavelength  streamwise  instabilities  necessarily  control  the  process; 
the  spanwise  instabilities,  paced  by  the  stage  of  development  of  the  supporting 
large  scale  structures  in  streamwise  planes,  then  follow,  and  just  respond  to, 
that  development.  Thus,  the  study  of  large  scale  structure  dynamics  in  a  two- 
dimensional  approximation  is  indicated  again.  In  this  approximation  we  then 
approach  the  specific  problems  posed  in  the  introduction. 

*The  measurements  of  Wygnanski  et  al  (1979a)  indicate  much  larger 
spanwise  correlation  lengths  for  a  velocity  ratio  A  =  0.4. 
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In  order  to  assess  what  constitutes  a  realistic  description  of  the  history 
and  structure  of  nonlinear  wave  motions  associated  with  individual  finite 
amplitude  Instabilities  we  consider  the  nonlinear  time  development  of 
homogeneous  Invlscld  Incompressible  two  dimensional  unstable  shear  flows. 
With  reference  to  cartesian  coordinates  (x,  y),  we  assume  that  the  Initial 
state  of  motion  results  from  the  superposition  of  a  single,  small  amplitude, 
excited  eigenmode,  periodic  In  x  and  t  (with  wavenumber  a  and  frequency  8) 
upon  a  mean  (x-averaged)  parallel  velocity  field;  that  Is,  at  t  ®  01 


V.(x,  0)  =  CU(y,0)  +  u  (x,  y,  0),  v(x,y,0)]  (1] 

We  describe  the  time  evolution  of  the  flow  In  terms  of  the  stream  function 
5l(x,y,t)  which  must  satisfy  the  vorticity  conservation  equation 

(v^)t  +  $y(V*$)x  -  $x  (V*$)y  *  0  (2) 

subject  to  the  initial  conditions  (1)  and  to  the  boundary  conditions 

y  -  ±  “  %  -  0  (3-a) 

if  the  flow  is  unbounded,  or  thejboundary  conditions 


y  =  ±  y, 


*x  =  0 


(3-b) 


if  the  flow  is  bounded  by  rigid  walls  at  y  =  ±  y 

w*  i 

We  seek  solutions  of  equation  (2)  under  the  assumption  that,  for  all  times 
t  >  0,  the  stream  function  $  may  be  resolved  into  a  mean  flow  contribution 
^(y.t)  and  a  disturbance  contribution  eij/(x,y,t),  the  latter  possessing 
small  but  finite  amplitude  characterized  by  the  parameter  e«l,  as  well 
as  x-periodicity  with  wavelength  (2tt/cx);  thus 

iHx.y.t)  =  ¥(y,t),  +  e^(x,y,t)  /a' 
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Substitution  of  (4)  into  (2)  yields 


(V2Y  +E.V2i())t  +  e  vt  -  i|»x  (V27)  ] 


+  e2  tyyTVV)x  -  ipx  (v24»)y3  -  0  | 


and  upon  x-averaging 

(7*',¥)t  +  ,e*<*  (7*  i|»>x  -  \|>X(V2  \p)  >  =  0 


which  is  the  equation  governing  the  time  evolution  of  the  mean  flow.  The 
equation  describing  the  evolution  of  the  periodic  disturbance  follows  readily 
upon  subtraction  of  (6)  from  (5),  viz. 

(V2^)t  +  v y(V2i|>)x  -  ^(7*7)  \ 


+e[4»y (V2xp)x  -  ^X(V2^ )y  -  <ty(72i|»)x  x(724»)y>3  s  0 

Equation  (6)  indicates  that,  if  the  disturbance  stream  function  remains  ! 
of  order  unity  under  differentiation  with  respect  to  the  space  coordinates, 
the  mean  flow  vorticity  (v2f)  varies  on  a  slow  time  scale  t=e2t.  Under  these 
conditions  the  solution  of  (7)  may  be  sought  by  the  method  of  multiple  scales 
(e.g.,  Nayfeh  1973),  and  a  uniform  expansion  for  the  disturbance  stream  func¬ 


tion  in  powers  of  e  may  be  obtained  according  to 

3  n 

e¥(x, y»t.)  =  s  en  i|>n(x,y,t,r)  ,  (8) 

n  =  1  / 

If, in  line  with  our  assumption»a  single  eigenmode  is  excited  at  t  =  0,  the 

leading  term  of  the  solution  (8)  must  have  the  form 


V>1(x,y,t,T)  =  1  /2 [A-j (x)  ^(y.x)  eQ  +  c.c.j 


where 


■» . ..v.’r:..  . : - * — * . . . 


'OICM 
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Substitution  into  (7)  then  shows  that  the  function  c|(y,x)  must  satisfy  the 
Rayleigh  equation 

L^)  a  (U a  >8)  (Ctyy-w2*])  -  o»Uyy  ^  =  0,  (10) 

with  the  homogeneous  boundary  conditions 

y  +  ±  00  5,  0;  (11) 

1  I 

but  that  the  function  A-j(t)  remains  undetermined  at  this  level  of  approxima¬ 
tion.  The  equation  governing  the  presumed  slow  time  evolution  of  the  mean 
flow  to  this  leading  order  is  obtained  readily  from  (6)  and  (9),  viz 

UT  =  -1/2  alAjJf/  exp(2  |  dt)  Im^C-yy)  (12-a) 

k 

where  denotes  the  complex  conjugate  of  c-j  and  B.j  *  Im(p).  If  (10)  is 

k  k 

used  to  express  £-jyy  in  terms  of  U  and  £j,  (12-a)  may  be  recast  vin  the  form 

£ 

Ut  “  f*  l*il*/exp(2/  B,*)  8,  (Uc.  -  (12-b) 

which  clearly  is  of  parabolic  type,  with  space  and  time  dependent  effective 
diffusivity,  as  long  as  8.  >  0.  Successive  differentiations  of  (12-a,b) 

with  respect  to  y  readily  yield  expressions  for  the  slow  time  variations  of 
the  slope  and  curvature  of  the  mean  velocity  profile  U(y,r),  viz. 

V  U.lAj^ptZ  /  »,*>  l">(C1C|yy)y  . 


2  /lAjixpU  /  B  dt)  B,  |Uct 
;  0  1  * 


+  V5i 


2a2lUa  -8|“2  U  U  U 

y  t 


(13) 
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l 

1,1,  t 

Uyyx  "  "  ?aiAil  exp  &  f0  *\  dt)  Im  ^lclyy}yy  * 

i  t  _ 

a2  |A,|i2exp  (2  /  p,  dt)  p,  jUo  -  B|  * 

T~  1  o  *  1 

1 


tUyy|Cl  |2yy  +  2  Uyyy  ICjly  +  Uyyyy  |C*|  ‘|  )  (W) 

-2a2  jUot  -  pf*  I(UUy)y  UT  +  2  U  UyUyT] 

Clearly  equation  O0)  and  (12)  merely  reproduce  the  model  of  Meksyn  and 
Stuart  (1J51)  for  the  study  of  finite  disturbances.  The  derivation  by  the 
multiple  scales  method  clarifies  the  crucial  underlying  assumption  of  slow 
variation  of  the  mean  flow  as  well  as  provides  the  framework  for  a  systematic 
evaluation  of  the  function  A-j(t)  (see  below)  without  appeal  to  the  heuristic 
energy  considerations  adopted  in  the  original  model.  The  equations  (13)  and 
(14)  have  no  direct  role  in  the  determination  of  the  flow  according  to  the 
model.  They  are  presented  here  because  experience  with  solutions  of  the 
Rayleigh  equation  (10)  Indicates  that  the  eigenfunction  t.j  and  the  associated 
eigenvalue  3  are  quite  sensitive  to  the  curvature  Uyy  of  the  mean  velocity 
profile.  An  examination  of  (14)  then  provides  the  key  to  the  a  posteriori 
assessment  of  the  renditions  where  slow  time  variation  of  the  mean  floisand 
attendant  description  of  the  disturbance  evolution  according  to  the  multiple 
scales  method, may  realistically  be  assumed.  Toward  that  objective,  we  need  to 

i 

determine  the  equations  governing  the  second  order  o(e2)  and  third  order  o(es) 
contributions  to  the  solution  (8);  this  is  done  below. 

Elaboration  of  (7)  in  the  context  of  the  solution  (8)  shows  that  the  contribu¬ 
tion  i|>2  to  the  disturbance  stream  function  must  have  the  form 


i 

i 


•;-3 


j 


J 

3 

J 


1 
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V2(x,y,t,T)  •  {  £A| (t )  Cg (y ,t )  e210  +  c.c.] 
the  function  Cg  being  the  solution  of  the  equation 

L2(c2)  =  (to  -  e)  (?2yy  -  4a2  c2)  -  allyy  c2  - 


a  ,  . 

"2"  1  clyyy  "  ?ly  ^lyy' 

subject  to  the  boundary  conditions  ' 

y  *►  ±  «  c  j  0  1 

However,  the  function  A^t)  still  remains  undetermined.  As  usual  In  weakly  j 
nonlinear  stability  theories,  an  equatlcii  for  Aj(t)  is  obtained  upon  consider-  I 
ation  of  the  contribution  o(e»)  to  the  solution  (8). Elaboration  of  (7)  to  d{es)  \ 
shows  that  thef  Inhomogeneous  part  iof  the  eq’tation  for  Includes,  among  i 

others,  terms  proportional  to  exp(10),  which  are  contributed  by  the  derivative  I 
of  with  respect  to  the  slow  time  variable  t.  Secular  terms  In  the  particular 
solution  for  ij»3  ,  and  the  attendant  nonuniformity  of  the  expansion  (8),  can  1 
then  be  avoided  only  by  seeking  the  particular  solution  in  the  form 


(16) 


(17) 


^(x.y.t.r)  -  j  [  Aj(t)  | Aj(t ) 1 2c3(y ,t)  e10  +  c.c.] 
the  function  c3  being  the  solution  of  the  equation 

i 

S(c3>  =  0J  «  -  3)  (Z3yy  -  azc3)  -  aUyy;3  = 


(18) 


1(Aj  lAil2)"1(A1g1)T+  exp!  i (9-  e*)  ]  g2 

(19) 

where 

g  i  (y»T) s  (clyy  -o2C2) 

(20-a) 

92 (y .T )  ?2y  C*wy)  -  2a  C,‘  (t2yy-3a!c2)- 

aCf(C2yy  -  3a2C2)y ]/4 

(20-b) 

-13- 
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Clearly  integrabillty  of  (19)  requires  that  J 

«o 

f  g1<t  -  1  AjlAjl^xp  [1(6-  6*)  ]  g2ftj  dy  »  0 

-00 

where  ^(y.x)  denotes  the  solution  of  the  equation  adjoint  to  the  Rayleigh 
equation,  viz 

(U  a  -  8)  &m  -  a^j)  +  2  aUy  *  0  (22) 

Accordingly  A^x)  must  satisfy  the  ordinary  differential  equation 

hj  AJ  +  h2  Aj  +  h^  A1|A1|2  =  0  (23) 

where 

00  ? 

h,  (x)  -f  g^y.x)  Cj(y»x)  dy  (24-a)  j 

A  -00  3 

00  ^ 

h2  (x)  =/  [(3/8x)  g^y.x)]  ^(y,x)  dy  (24-b)  | 

—oo  4 

at  4 

h3  (x)  =  -  1  exp  [i  (6  -  8*)]/  g2(y»x)  &  (y.x)dy  (24-c)  | 

^  *.00  -jj 

Equation  (23)  embodies  the  extension  of  the  Meksyn  and  Stuart  model  yielded 
by  multiple  scales  considerations.  Interestingly  enough.  If  the  x  dependence  i 

j 

of  the  eigenfunction  is  Ignored,  i.e. ,  if  glT  *  h2  s  0  in  (24-b),  equation  (23)  j 

reduces  to  the  classical  result  of  Stuart  (1960);  thus  both  the  Meksyn  and  Stuart  j 

J 

and  the  Stuart  models  arise  as  particular  cases  in  the  present  context.  ,  i 


The  derivative  gjT  in  (24-b)  and  the  coefficient  h2  in  (23)  are  retained  here;  they| 
are  evaluated  ty  noting  that,  since  the  eigenfunction  r,  satisfies  the  Rayleigh  j 
equation  (10),  the  derivative  must  satisfy  the  equation  j 
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readily  derivable  from  (10)  upon  differentiation  with  respect  to  the 
slow  time  t. 

Equation  (25)  Is  also  subject  to  an  <ntegrab1 flty  condition  which  yields  an 
ordinary  differential  equation  for  the  t  dependent  eigenvalue  6(t),  namely 

Pr“tb4(t)/h1(x)]  or  j  &(t)  -6(0  W  Oyhj)  dr  (26) 

o 

where  \ 

VT>  •-«£  <ly  (27) 

With  this  condition,  which  merely  reflects  the  assertion  between  the 
T-dependent  eigenvalues  6  an».  elgensolutlons  Cj  of  equation  (10),  the  formu¬ 
lation  of  the  extended  Meksyn  and  Stuart  model  is'  completed.  Its  applica¬ 
tion  to  determine  the  evolution  of  an  Invlscldly  unstable  free  shear  flow 
subject  to  finite  amplitude  Instabilities  involves  a  stepwise  advancement  of 
the  solution  in  time.  Given  the  initial  mean  velocity  field  U(y,  0)and  the 
wavenumber  a  as  well  as  the  initial  amplitude  Aj(0)  of  the  disturbance,  the 

j 

unit  step  In  that  process,  between  times  Oand  At,  Involves  sequentially: 

1)  the  integration  of  equations  (10),  (22)  and  (16)  to  determine,  respectively, 
the  eigenvalue  8(0)  together  with  the  function  c-(y,0)  and  the  functions 

i 

^'(y.o),  C2(y.O);  1i)  the  evaluation  of  UT(y,0),  and  U  (y  0),  8'(0)  from 
(12),  (14)  and  (26),  respectively,  and  the  consequent  determination  of 
ClT(y»0)  by  integration  of  (25);  Hi)  the  integration  of  (23)  and  (12-b)  to  ob¬ 
tain  Aj(At)  and  U(y,Ax),  whereupon  the  next  step  can  be  started.  Numerical  pro¬ 
cedures  for  these  calculations  are  presented  in  Section  3,  and  an  assessment  of 
the  attendant  predictions  for  the  evolvina  structure  of  the  wave  motion  associated 
with  a  finite  amplitude  instability  is  set  forth  in  Section  4. 
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Whereas  measurements*  even  If  detailed  and  spectrally  resolved*  do  not 
usually  provide  much  Information  about  those  features  of  nonlinear  Instability 
waves  —  viz. ,  transversal  variation  In  phase  and  associated  velocity  as 
well  as  vortlclty  fields  —  which  are  most  Important  for  the  realistic  evalu¬ 
ation  and  refinement  of  a  theoretical  model*  the  assessment  of  Section  IV 
Includes  a  three  way  comparison  among  available  experimental  data*  predictions 
of  the  (weakly  nonlinear)  extended  Meksvn  and  Stuart  model,  and  predictions 
of  a  strongly  nonlinear  model.  The  choice  of  the  nonlinear  model  Is  tailored 
to  both  lines  of  Inquiry  reported  1  r.  the  paper,  namely:  ythe  identification 
of  the  cyclic  subharmonic  generation  mechanism  In  mixing  layers  and  the  descrip¬ 
tion  of  the  detailed  structure  of  the  nonlinear  wave  motions  associated  there¬ 
with.*  As  stated  above*  we  seek  to  demonstrate  that  the  dynamics  of  two- 
dimensional  inviscidly  unstable  transitional  and  turbulent  shear  flows  are 
dominated  by  the  same  deterministic,  self -preserving,  repetitive  process  of 
subharmonic  generation  wherein  the  nonlinear  evolution  of  an  individual 
finite-amplitude,  two  dimensional, streamwise-propagating  instability  intrin¬ 
sically  provides  mean  flow  growth  as  well  as  cyclic  excitation  of,  and  sub¬ 
sequent  energy  transfer  to, its  own  1/2  subharmonic.  In  this  view  .a  re¬ 
stricted  class  of  nonlinear  Instability  processes  is  singled  out  as  govern¬ 
ing  the  development  of  turbulent  mixing  lexers.  Specifically,  the  initial 

t 

formation  of  vorticity  concentrations  In  the  transition  region  is  associated 
with  the  self-limiting  nonlinear  growth  of  a  single  dominant  inviscid  insta¬ 
bility  —  the  fundamental  --  whose  selection  is  governed  by  a  convolution  of 
the  disturbance  environment  with  the  receptivity  (Tam  1978)  and  the  frequency/ 

*  The  failure  of  the  weakly  non! inear  model  in  describing  the  structure 
of  Individual  instabilities  (see  Section  4)  preempts  its  extension  to 
the  subharmonic  generation  problem. 
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amplification  rate  properties  (Mlohalke  1964)  of  the  Initial  mixing  layer. 

The  evolution  of  the  fundamental  Is  strongly  Influenced  by  the  concurrent 
distortion  of  the  mean  flow,  but  is  largely  unaffected  by  the  concurrent 
excitation  of  harmonics  (Zabusky  and  Deem  1971).  Upon  saturation  the 
fundamental  undergoes  a  self-induced  modulation  (the  vortex  nutation  ob¬ 
served  by  Zabusky  and  Deem  1971,  Patialk  et  al  1976,  Miura  and  stt  o  1978), 
which  affects  the  amplitude  and  the  structure  of  the  wave  motion,  as  u-ell 
as  the  distortion  of  the  mean  flow.  As  we  shall  demonstrate  In  section  5,  the 
frequency  of  this  modulation  Is  nearly  1/2  that  of  the  saturated  fundamental, 
while  the  concurrent  mean  flow  distortion  is  such  thfct  the  subharmonic  fre¬ 
quency  closely  approaches  those  of  the  maximum  amplification  rate  spatial  and 
temporal  Instabilities  for  the  prevailing  mean  flow.  Since  the  phase  veloci¬ 
ties  of  these  instabilities  and  of  the  fundamental  are  then  essentially  the 
same.  1/2  subharmonic  excitation,  with  well  defined  initial  amplitude  as  well 
as  structure  of  the  attendant  wave  motion,  ensues  from  the  modulation  of  the 
fundamental.  The  intrinsically  excited  subharmonic  then  amplifies,,  draining 
energy  from  the  fundamental  as  well  as  from  further  distortion  of  the  mean 
flow.  However,  its  growth  is  again  self-limiting  and  such  that,  at  satura¬ 
tion,  th?  original  fundamental  becomes  nearly  spent,  while  the  structures  of 
the  mean  flow  and  of  the  subharmonic  wave  motion  reproduce,  except  for  a  doub¬ 
ling  of  all  length  scales,  the  structures  that  prevailed  at  the  earlier  time 
of  fundamental  saturation.  A  self-preserving  process  of  cyclic  nonlinear 
subharmonic  generation  and  mean  flow  development  follows  necessarily. 

Whereas  only  two  waves  (a  fundamental  and  its  second  harmonic  interacting  with 
each  other  as  well  as  with  the  mean  flow)  are  active  at  any  stage  of  the  non_ 
linsar  wave  cycle  postulated  above,  our  nonlinear  model  seeks  solutions  of  the 
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equations  (6)  and  (7)  under  the  assumption  that  the  perturbation  stream 
function  tp  may  be  described  at  all  times  by  a  Fourier  series  In  x  truncated 
at  the  second  term,  viz 

2 

Mx,y,t)  -  imz  <>  (y,t)e1nax  +  c.c.]  ^ 

n-1  n 

With  this  assumption  the  modal  structure  of  the  evolving  flow  Is  clearly 
manifested  at  all  times,  and  the  complexity  as  well  as  cost  of  the  numerical 
Integrations  are  considerably  reduced.  In  Its  support  we  note  that  the  numer¬ 
ical  results  of  Mlura  and  Sato  (1978),  who  adopted  the  Fourier  series  approxi¬ 
mation,  and  those  of  Zabusky  and  Deem  (1971),  who  directly  Integrated  the 
time-dependent  incompressible  Navi er- Stokes  equations,  exhibit  qualita¬ 
tively  Identical  saturation  and  post-saturation  behaviors  for  finite  amplitude 
Instabilities  in  free  shear  flows.  Thus,  the  departure  of  finite  amplitude 
instabilities  from  an  Initially  sinusoidal  spatial  shape  (Patnalk  et  al  1976) 
seemingly  plays  a  secondary  role  in  their  development.  However,  further 
scrutiny  may  be  In  order  as  to  the  influence  of  that  departure  on  the  onset 
and  evolution  of  small  scale  secondary  instabilities  --  e.g.,  those  respon¬ 
sible  for  the  skewness  and  apparent  branching  of  the  mixing  layer  struc¬ 
tures  —  which  are  sensitive  to  the  details  of  the  flow  pattern  relative  to 
those  very  struct  ves. 

Substitution  of  the  assumed  form  of  solution  (28)  into  (6)  and  (7)  readily 
yields  the  system  of  coupled  nonlinear  partial  differential  equations  which 
govern  the  concurrent  time  evolution  of  the  mean  flow  and  the  selected  two 
dimensional  disturbance  modes,  namely 

“t  ■',Tt*l*!»y+i*lyy  +  2(*2  +2y>  -  **2  *2yy)l  j  <29-»> 
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t\mt)  -  UoOtoiyy  -  aUyy  • 

-«2  [2<l»2  (♦tyy  "  “  *l(*2yy  '  4  “^y  +  / 

♦2y  4lyy  -  aa*l>  ’  2  ♦ty^Zyy  “  4  «J W  / 


mmt)  -  2Ual  .  (♦gyy  -  4  a24>2)  +  2aUyy  ^  •  i 
-5“f  l*l(*lyy  ■  a*«>l>y  “  ♦ly^lyy  ‘  'V1  I 


(29-c) 


to  be  Integrated  numerically  subject  to  prescribed  Initial  conditions 
M(y .0) ,  4>1(y *0) »  $2(y,0)  and  to  J.he  boundary  conditions 
y  ±  «*  U(y,t)  ♦  U(y,0) 

♦i(y»t)  -*•  Q  <t>2(ytt)  *■  o 


Procedures  for  this  numerical  Integration  are  presented  In  the  ^following 
section,  and  attendant  results,  relevant  to  our  Inquiries,  are  reported 
and  discussed  In  Sections  4  and  5. 
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3.  NUMERICAL  PROCEDURES 

As  pointed  out  In  Section  2,  the  unit  problem  of  determining  the  evolution 
of  the  motion  associated  with  a  nonlinear  instability  over  a  time  Interval  At 
according  to  the  extended  Meksyn  and  Stuart  model  Involves  the  sequential  Inte¬ 
gration  of  equations  (10),  (22),  (16),  (25),  (23)  and  (12-b).  The  Integration 
of  equation  (10)  to  determine  Cj(y,t)  poses  a  classical  eigenvalue  problem, 
which  Is  solved  by  Iteration.  Given  a  trial  eigenvalue  B(t),  two  linearly 
Independent  solutions  and  satisfying  the  asymptotic  behavior 

Cj^(t)  exp  (-<xy)  for  y  -*■  *,  are  determined  by  numerical  Integration 
using  a  variable  order  Adams  predictor-corrector  method  to  march  from  large 
y  toward  y  »  O',  where  the  condition  [Re(C]^)/Re(cj^)  1 j  «  [Imfc^)/Im(Cj^)l 
is  to  be  satisfied.  The  value  of  8  Is  Incremented  using  a  Newton-Raphson 
scheme,  and  the  numerical  Integrations  are  repeated,  until  the  condition  on 
y  a  0  is  satisfied.  The  final  solutions  are  linearly  combined  so  that 

Re  Icly(0.T)]  c  Im  [^(O.t)]  «  0,  and  then  normalized  so  that  Re  [^(O,  t)1  a  1 
Exactly  the  same  procedure  is  then  followed  to  determine  {^(yyt)  by  numerical 
Integration  of  (22);  however,  no  Iteration  in  B  Is  required  because  (22) 
is  adjoint  to  (10).  The  determination  of  c2(y,t)  by  integration  of  the 
nonhomogeneous  equation  (16)  subject  to  the  homogeneous  boundary  conditions 
(17)  is  uniqueiand  reducible  to  the  evaluation  of  the  particular  integral, 
because  no  complementary  functions  of  (16)  exist  which  satisfy  (17).  Two 
particular  integrals  e2^  and  C2^»  satisfyin9  the  asymptotic  behavior 
5  C2^  exp(-2ay)  for  y -r  are  first  determined  by  numerical  Inte¬ 

gration  using  the  aforenoted  Adams  method  to  march  from  large  y  to  y  =  0;  the 
final  solution  is  then  obtained  by  the  linear  combination  c2  =  + 

(1  -  D2)  ?2^,  which  satisfies  the  conditions  Re[c2y(0,T)l  ■  lm[c2(0,x)l  *  0 


\ 
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Dj  ■  -A'‘[(C|1)*-c|2)*)  Re(;«>)+ 1  (C^)*  -52<2)*)  Im  (c2<8>)]  .  „ 

with  .  '  ,  • 

A  -ReKcl11  -  c|2))  (C^1**  -  ?sj2)*))  y  .  „ 

The  Integration  of  (25)*  subject  to  the  boundary  conditions  O.for  y  -►  » 

®nd  C^T  *  0  for  y  *  0,  is  reduced  to  the  evaluation  of  quadratures  by  the 
substitution 

y 

clT  (y.T)  »  e^y.x)  /  n,  (y,t)  dy 
1  0 

which  transforms  (25)  Into  the  equation 
n1y  +  2cl"X  ?ly  nl  *  f^»T^ 

with 

f(y,T)  *<Ua  -er*  (aUyyTCj  -  (UT  a-B')  (c^  -  a1^)) 

The  solution  is 

n1(y»t)  =  /  Cj  f  dy 

00 

Upon  this  determination  of  clT(y,x)  the  amplitude  function  AjU  +  Ax)  Is 
calculated  from  (23)  and,  finally,  the  mean  velocity  U(y,x  +  Ax)  is  obtained 
by  numerical  solution  of  the  nonlinear  parabolic  equation  (12-b).  In  that 
connection  a  standard,  stable,  implicit,  second  order  accurate  scheme, 
where  time  and  space  derivatives  are  approximated  by  centered  differences, 
is  used  to  convert  (12-b)  into  a  system  of  simultaneous  linear  equations 
for  the  unknowns  w.n+1  =  (Ua0*1  -  uj)  a  [Ua  (x  +  Ax)  -  UL(t)] ,  where 
j  55  1,2...  N  denotes  the  general  y  -  mesh  point.  When  the  boundary  conditions 
are  taken  into  account,  the  system  takes  tridiagonal  form  amenable  to  rapid 
direct  solution.  That  completes  the  advancement  of  the  weakly  nonlinear, 
extended  Meksyn  and  Stuart  solution  by  a  unit  time  step. 
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V 


Passing  to  the  nonlinear  model  described  at  the  end  of  Section  2,  the  unit 
problem  of  determining  the  interactive  evolution  over  a  time  Interval  At  of 
the  mean  flow  and  the  wave  motions  associated  with  a  nonlinear  Instability 
as  well  as  Its  bound  second  harmonic  Involves  the  simultaneous  Integration 
rtf  equations  (29-a,b,c),  In  that  connection,  an  impllcily,  second  order 
eccurate  scheme,  where  the  time  and  space  derivatives  are  approximated  by 
centered  differences,  is  used  to  convert  (29-a,b,c)  into  a  system  of  simul¬ 
taneous  linear  equations 


fln,n+l  j  nn*n+l  *  rn*n+1  4.  r>n*n+l  4.  cfiAlftl  -  *n 

Aj  *j-2  *  Bj  *J-1  +  Cj*j  +Dj*W  E/jV2-fo 


**  / 


(31) 


where  j  =  1, 2. .  .N,*!)+1  denotes  the  vector  with  components  (u!?+1-  U*?» 

J  J  J 

^lj^lj’  ”^2j ^  t0  be  determined  at  the  general  y-mesh  point  j,  and 


i 

■W 

■A 


the  coefficients  A?  through  E?  are  5x5  matrices.  When  the  boundary  conditions 

J  J 

(viz  ^0=  are  taken  into  account,  the  system  of  equations  (31) 

takes  block  pentadiagonal  form  and  thus,  becomes  amenable  to  direct  solution. 
The  algorithm  remains  quite  stable  when  the  time  step  At  is  chosen  consistent 
with  the  condition 

At 


max 


trjt|U|j +e|<f>1ylj  +e|4*yC3+  e«C|4>  In  +21*  i"}]  1 

0  1  J  2  J 


where 


ro 


eatl^iyy  +zi*2yy~  4ot2<Mj] 

lUyy'j  l*i www"  «**iwU  +  H,MMw***! 


lyyy 


ly'j 


2yyy 


V?  ■ 


approximately  describes  the  ratio  between  the  x-  and  y-derivatives  of  the 
vorticity  at  the  jth  mesh  point  and  time  t. 
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4.  EVOLVING  STRUCTURE  OF  NONLINEAR  UNSTABLE  HAVE  MOTIONS 
We  consider  the  evolution  of  a  temporally  developing  shear  layer.  The  Initial 
state  of  motion  is  characterized.  In  dimensionless  variables,  by  the  mean 
velocity  U(y,0)  =  tanh  2y,  upon  which  the  velocity  and  vorticlty  fields  atten¬ 
dant  to  the  nearly  most  unstable  eigenmode  a »  0.8892,  6(0)  *  0  +  0.3784  1  are 
superposed  with  s=  0.10  and  Aj(0)  *  0.10.;  Two  rigid  plates  located  at  y  *  ±3.875 
are  assumed  to  bound  the  flow.and  solutions  are  determined  subject  to  the  atten- 
dant  boundary  conditions. 


Analysis  of  the  flow  development  according  to  the  extended,  weakly  nonlinear, 

Meksyn  and  Stuart  model  yields  the  suspicious  result  that,as  time  progresses,  6^ 

decreases  but  A,  increases  (Figure  1)  to  such  an  extent!  that  the  amplitude 

t 

a(t)  =  cA,(t)  exp  (/  p^dt)  of  the  perturbation  stream  function  on  the  axis  y  a  0 
1  o  2 

distinctly  exceeds  the  linear  prediction  eAj(O)  exp  (0j(d)t)  when  a(t)  ^(°) 

(Figure  2).  Whereas  this  prediction  differs  from  that  of  the  nonlinear  model 
(Figure  2),  as  well  as  from  the  trends  suggested  by  measurements  in  transi¬ 
tional  mixing  layers  forced  at  the  most  unstable  frequency  (Figure  3),  a 
failure  of  the  Meksyn  and  Stuart  approach  Is  indicated.  A  little  diagnostics 
readily  locates  the  source  of  the  difficulty  in  the  assumption  of  slow  time 
variation:  even  at  the  initial  time  the  derivative  (clyy  "  a2Cj)T  of  / 
order  unity,  much  larger  than  ^ (Figure  4). 


The  trend  of  (clyy  -  a2£|.)T  versus  y  shown  in  Figure  4  further  suggests  that 

the  problem  stems  from  an  inadequate  description  of  the  vorticity  transports 

in  the  cat's  eye  of  che  evolving  wave  around  the  critical  layer  y  =  0.  In 

t 

that  region,  having  time  dependent  lateral  extent  6(t)  exp(/  0-jdt)/ 

Uv(Q,t)]1//2  =  [a(t)/U  (0,t)]  ly/2,  the  mean  velocity  U  as  well  as  the  curvature 
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\ 


of  the  mean  velocity  profile  U  are  o(6);  however,  by  virtue  of  (12-t)  and 

«/  J 

(14),  the  time  rates  of  change  of  those  quantities  when  aUy(0,x)6(t)>Bi(T) 
become  respectively,^"^  -  o(6s)  and  &-j"luyyt  =  o(6).  Under  those  conditions 
the  curvature  Uyy,  which  has  a  profound  Influence  on  the  shape  of  the  eigen¬ 
function  exhibits  in  the  critical  layer  a  normalized  rate  of  change 
^iUyy^_1  ^yyt  =  comparable  to  that  of  the  disturbance  vorticity,  viz 
[BiCCjyy  -  a2^)]  1  (?iyy  -  «2?1)t  *  6  (1).  As  a  consequence,  the  assumption 
that  the  structures  of  the  mean  flow  and  of  the  vorticity  field  associated  with 
the  instability  vary  on  a  slow  time  scale  collapses,  and  the  weakly  nonlinear 
model  fails  to  describe  the  flow  realistically.  The  extent  of  such  failure 
is  vividly  demonstrated  in  Figure  5,  where  the  structures  of  the  wave  motion 
predicted  by  the  extended  Meksyn  and  Stuart  model  and  by  the  nonlinear  model 
for  t  =  9.82»a(t)  =  .285, 3(x)  =  .1767  are  compared  with  each  jther  as  well  as 
with  the  structure  prevailing  at  t  =  0.  Even  greater  divergences  arise  in  a 
comparison  between  the  growth  rates  and  the  structures  of  the  wave  motion 

Dredicted  bv  the  two  models  for  the  second  harmonic,  i.e..  the  normalized 
S2(y.T)/RelC2(0,x)J  and  its  derivatives  vis  a  vis  the  normalized 

<i>2(y.t)/Re[<i>2(0,t)]  and  its  derivatives.  Thus,  reliance  on  the  nonlinear 
model  is  indicated,  and  justified  below,  for  the  problem  on  hand. 

) 

The  evolving  structure  of  the  wave  motion  which  the  nonlinear  model  predicts 
for  the  amplifying  instability  up  to  the  time  of  its  saturation  (t  $  13.5 
in  Figures  2  and  6)  compares  favorably  with  experiments  on  at  least  two 
counts.  First  is  the  qualitative  !  agreement  of  the  relative  growth  rate 
trends  predicted  by  the  model  for  the  scaled  energies 

Etn(na,t)  =/W  e2(|$ny|2+  n2a2|<f>J2)  dy 
-yw  ? 
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of  the  dominant  fundamental  (n^l)  and  its  second  harmonic  (n=2)  In  the  time 
developing  shear  layer  Identified  at  the  beginning  of  this  sect1on(F1gure  6) 
vis  a  vis  the  relative  trends  of  the  energies 

r.m.s.  (n8,x,y)  dy 


W-M)  -/*  “  ,  . 

%  2^<xo)  UI 


measured  by  Mlksad  (1972)  In  a  space  developing  transitional  mixing  layer  forced 
at  the  same  scaled  frequency  (Figure  3):  In  both  cases  the  energy  of  the  har¬ 
monic  grows  at  a  rate  approximately  1.5  times  that  of  the  fundamental,  a  be¬ 
havior  certainly  not  predicted  by  weakly  nonlinear  theories.  We  postpone  analy¬ 
sis  of  this  behavior,  anticipated  qualitatively  by  Robinson  (1974)  and  dis¬ 
cussed  further  below,  to  elaborate  the  second  count,  viz.  the  quantitative 
agreement  between  ‘•.he  magnitudes  of  the  energies  Etl  and  Eul  which  theory  and 
experiments,  respectively,  associate  with  the  fundamental  at  saturation.  Upon 
reference  to  the  scales  of  the  equivalent  temporally  developing  mixing  layer, 
we  set  8m(xQ)  *  1/4  and  Uj  =  2  in  the  definition  of  Eul*  We  then  proceed  to 
calculate  the  ratio  (Et i/Eui)  consistent  with  these  values  and  with  the  theor¬ 
etically  predicted  structure  of  the  wave  motion  at  saturation.  The  result  is 
(Eti/Eui)  ®  3. 5,  a  value  nicely  recovered  upon  comparison  of  Figures  3  and  6. 


From  figures  3  and  6  we  also  see  that  the  energies  E^  and  Eu2  associated 
with  the  second  harmonic  compare  favorably  during  the  early  phases  of  devel¬ 
opment,  but  depart  from  each  other  as  the  fundamental  approaches  saturation. 

A  spurious  double  hump  arises  in  the  theoretical  predictions  at  that  stage. 

We  discount  the  significance  of  the  hump  feature  In  the  assessment  of  the 
nonlinear  model  on  three  counts,  namely:  i)  the  energy  associated  with  the 
harmonic  remains  small  compared  to  that  of  the  fundamental  and  exhibits 

*1 

I 
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growth/dtcay  trends  apparently  controlled  by  the  evolution  of  the  fundamental 
Interacting  with  the  mean  flow  (Figure  3);  11)  the  qualitative  agreement  be¬ 
tween  the  results  of  Zabusky  and  Deem  (1971)  and  those  of  Miura  and  Sato 
(1978)  confirms  the  secondary  impact  of  the  harmonic  upon  the  saturation  and 
post-saturation  evolution  of  the  fundamental,  i.e. ,  on  the  processes  which 
underlie  the  subharmonic  generation  mechanism  proposed  and  addressed  In  this 
paper;  ill)  the  double  hump  feature,  which  repeatedly  arises  in  calculations 
using  a  truncated  Fourier  series  representation  of  nonlinear  wave  motions 
approaching  saturation  (see  section  5),  is  eliminated  in  unsteady  twodimen¬ 
sional  flow  simulations,  which  do  not  introduce  the  truncation  approximation  m 
(see  Zabusky  and  Deem  1971).  Thus,  truncation  and  the  attendant  neglect  o' 
the  higher  harmonics  are  responsible  for  the  observed  divergence  between 
the  energies  anc!  Eu,>;  however,  the  approximation  does  not  obscure  the 

essential  aspects  of  nonlinear  instability  evolution  evidenced  by  the  calcu¬ 
li 

lations, and  indirectly  confirmed  by  the  experiments  as  discussed  below.  \ 

1 

The  prime  aspect  of  finite  amplitude  instability  dynamics  resides  in  the  pre-  ] 

vlously  noted, relatively  rapid, evolution  of  the  curvature  U  of  the  mean  i 

yy  -j 

velocity  profile  within  the  critical  layer  of  lateral  extent  fi(t).  We  have  s 

indicated  that  U  undergoes  rapid  change  when  and  if  the  nonlinear  Instability 

'J  J  1 

attains  a  small,  but  finite,  amplitude  a(t)  such  that  a  Uy(0,x)  6(t?  >  B^(t);  j 

conversely,  the  test  a  Uu(Q,x)  6(t)<Mx)  defines  the  conditions  under  which  i 

y  i 

weakly  nonlinear  models  may  be  employed  to  realistically  describe  the  structure  I 
of  the  unstable  wave  motion.  Beyond  that  threshold  equations  (13)  and  (14)  show  | 
that  the  overall  flow  evolution  becomes  characterized  by  a  time  rate  of  change  j 
of  mean  flow  vorticity  within  the  critical  layer  =  o(62)  *  oCa(t)3  j 

comparable  to  that  of  the  disturbance  amplitude  B^_1(da/dt)  =  a(t),  but  a  J 

time  rate  of  change  of  the  y-derivative  of  mean  flow  vorticity  B^-1  U  ^  *o(6)  | 
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larger  than  either  of  those  two  within  the  critical  leyer.  Clearly,  the 
total  vortlcity  Va$  In  that  layer  is  then  described  by 

V*iJ»(x,y,t)  *  Uy(0,0)  +  a(t)  ft^x.Y.t)  (32) 

where  fli  denotes  c.  function  which  Is  of  order  uniuy  together  with  its 
Y-derivatlve  with  respect  to  the  stretched  coordinate 


Y  -  [Uy(o,t)/a(t)]1/2  y  *  (y/6) 


(33) 


across  the  critical  lawyer  centered  at  y  *  0.  If  the  considered  instability, 
characterized  by 

eij>(x,y,t)  “  a(t)  [cosox  +  o(<5)] 

in  the  neighborhood  of  y  -  Q,  dominates  the  disturbance  field,  the  velocities 
within  the  critical  layer  Y  =  o( 1)  are  to  leading  order 

u  -  Uy(0,t)  H  (34-a) 

v  .»  a(t)as1nox  . . (34-b) 

Conservation  of  vorticity  then  requires  that  the  function  U.  defined  by  (32) 
satisfy  the  equation 

n1T  +  Yflie  +  sin  5  fliy  +  -  (1  +  i  niY3  =  0  (35) 

in  "boundary  layer"  coordinates  (£,Y,T)  defined  by 

(36) 


K  - 


y  *  (y/<sh 


dT  =  (ao/6)dt 


the  quantities 

w  (T)  =  (d  log  O  (T)/dT)  -  6,10*0)  Uy(0,T)l  *  “i [aV°-T)S1 ^  (37-«) 


V2(T)  =  -  [d  log  Uy(0,T)/dT]  j 


x37-bi 


being  of  comparable  magnitude  <  1  under  the  considered  conditions  (3i/«S)  <1. 
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Inspection  of  (35)  readily  suggests  that,  due  to  the  vortlclty  transport 
s1n£  OiY.  the  form  of  Its  solution,  periodic  In  £  over  the  wavelength  2n  of  a 
dominant  monochromatic  Input  disturbance.  Is  necessarily 

Q,(£,Y,T)  -  ft/0*  (Y,T)  +  1/21  I  fl (Y,T)fe1n5  +  c.c.l  (38) 

1  1  h*  1 

even  though  (Y,0)s  0  for  n  >2.  Thus,  hamonlcs  are  forced  by  the 

noted  transport  within  the  critical  layer.  Substitution  of  (38)  into  (35 ) * 
and  Inspection  of  the  attendant  equations  for  the  functions  '  (Y,T), 
readily  show  that  the  forced  vorticitles  a(T)  (Y,T),  n  >  2,  become  of 

order  a(T)  after  time  Intervals  T  =  o(l)  and,  by  that  process,  harmonic  pertur¬ 
bation  velocities  of  order  a^^(T)  are  generated  at  the  edges  of  the  critical 
layer.  Accordingly,  harmonics  observed  in  experiments  and  determined  in 
nonlinear  calculations  should  exhibit  amplification  rates  approximately  1.5 
times  as  large  as  that  of  the  dominant  fundamental.  This  Is, in  fact,  the 
trend  evidenced  by  the  results  of  Figures  3  and  6,  and,  more  generally,  by 
Miksad's  (1972)  observation  that  "the  measured  growth  rates  of  harmonic  modes 
range  from  1.15  to  1.73  times  that  of  the  fundamental"  in  transitional  mixing 
layers  monochromatic ally  forced  at  various  scaled  unstable  frequencies.  The 
scatter  of  the  observed  growth  rate  ratios  about  the  value  1.5  is  not  sur¬ 
prising  because  different  ihilstories  of  the  parameters  Pj(T)  and  p2(t)  pertain 
to  different  situations;  accordingly,  the  components  n]  y(Y,T)  of  the  solution 
(38)  vary  from  case  to  case  and  so  do  the  ratios  between  pairs  of  attendant 
forced  perturbation  velocities  and  eneVgies. 

The  important  role  of  nonlinear  vorticity  transports  within  the  critical  layer 

in  determining  the  behavior  of  finite  amplitude  instabilities  and  associated 

harmonics  was  recognized  by  Robinson  (1974)  who,  In  fact,  proposed  an  asymptoti 
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solution  to  a  reduced  form  of  equation  (35).  In  that  context.  Robinson 
reached  several  perplexing  conclusions,  which  contrast  with  experimental 
evidence  and,  thus,  warrant  specific  comment.  At  Issue,  are  the  conclusions 
that:  1)  the  temporal  growth  rate  of  Instabilities  with  wavenumber  close 
to  the  marginal  value  1  fon  a  tanhy  mean:  velocity  profile  should  be 
c.j  a  -  3(a-  1)  Instead  of  *-0.637(a-  1)  as  predicted  by  linear  theory 
and  Indicated  by  experiments  (Mlksad  1972),  at  least  In  the  early  phases 

4 

of  nonlinear  amplification;  11)  no  spatially  amplifying  solutions  should 
exist  In  the  nonlinear  range  even  though  data,  such  as  those  of  Figure  3, 

*  "3 

clearly  display  their  presence;  Hi)  the  harmonics,  forced  by  the  nonlinear 


redistribution  of  vortlclty  within  the  critical  layer,  should  amplify  at  a 
rate  exactly  1.5  times  as  large  as  that  of  the  dominant  fundamental .rather  than 
the  varying,  although  comparable,  rates  noted  and  rationalized  above.  Rejection 


of  these  conclusions  is  suggested  by  the  following  critical  review  of 
Robinson's  solution. 

A  change  of  independent  variables  from  (£,  Y»T)  to  (£,¥,  T),  where 

•»  v2 

¥  «-£  +  cos  s 

denotes  a  scaled  stream  function  In  the  critical  layer,  reduces  (35)  to 
the  equation 

°1.T+  ¥!!i?  +  “l  [tli  -  ni»  1  "  0 

In  turn  (39-a)  reduces  to  the  equation  (4.28)  of  Robinson,  namely 

Yn1C  +  wl  (ni  “  )=0 

under  the  assumptions/approximations  that 

a  =1,  Uy(0,T)  *  1,  y2(T)  =  n.T  (£,  T)  =0 
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whereupon  the  parameter  Vj  *  c^A"1^2  In  his  notation.  Robinson  determined 
an  asymptotic  solution  to  the  reduced  equation  (39-b)  for  «  1  subject  to 
the  boundary  conditions  that,  as  ?-*■«,  the  vortlclty  matches  a  time 
Invariant  mean  vortlclty  distribution  upon  which  Is  superposed  the  vortlclty 
of  a  linear  Instability  having  a  time  Invariant  wave  motion  as  well  as  a 
time  Invariant  amplification  rate  8^  Clearly  these  boundary  conditions  do 
not  describe  the  actual  state  of  affairs  for  even  the  weakly  nonlinear  Meksyn 
and  Stuart  model  predicts  that  neither  the  mean  flow  vortlclty  (Figure  7)  nor 
the  structure  (Figure  5)  and  amplification  rate  (Figures  1  and  2)  of  the  un¬ 
stable  wave  motion  remain  time  Invariant  outside  the  critical  layer.  In  fact, 
Robinson's  results  contradict  the  assumed  asymptotic  behavior  for  they  reveal 
a  velocity  jump  of  order  a‘3^2(T)  across  the  critical  layer,  a  jump  which  is 
assumed  to  be  absorbed  by  the  unspecified,  and  unfortunately  neglected, 

"(a)  spreading  or  contraction  of  the  basic  flow,  (b)  similar  spreading  or  con¬ 
traction  of  the  fundamental  disturbance,"  as  well  as  by  the  introduction  of 
all  higher-harminlcs  outside  the  critical  layer. 

The  inadequate  boundary  conditions  lead  to  an  asymptotic  solution  for  flj(S,Y,T) 
which,  according  to  Robinson's  equation  (4.48),  is  completely  defined  by  the 
instantaneous  value  of  the  parameter  ^(T),  independent  of  initial  conditions 
and/or  previous  evolution  of  the  disturbance.  Whereas  this  behavior  Is  incom¬ 
patible  with  the  indications  of  either  nonlinear  models  or  experiments  —  which 
unequivocally  show  an  influence  of  prior  history  upon  the  evolving  structure  of 
nonlinear  unstable  wave  motions  —  the  particular  nature  of  Robinson's  solution 
as  well  as  the  resolution  of  the  perplexing  conclusions  it  leads  to,  become 
apparent;  In  connection  with  the  conclusion  1)  we  note  that  the  redistribution 
of  vorticity  within  the  critical  layer  of  an  evolving  nonlinear  instability  must 
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be  determined  by  solving  equation  (35),  where  the  time  derivative  Is  re¬ 
tained,  and  not  the  Robinson's  equation  [(30-b)  In  this  paper],  where  Is 
neglected.  Retention  of  the  time  derivative  allows  the  structure  as  well  as 
the  amplification  rate  of  a  nonlinear  Instability  with  initially  small  ampli¬ 
tude  at  first  to  follow,  and  then  to  depart  from,  the  predictions  of  linear 
theory.  In  connection  with  the  conclusion  11)  we  note  that  the  strong  Inter¬ 
action  between,  and  the  concurrent  evolution  of,  the  mean  flow  and  the  wave 
motion  structures  require  that  a  realistic  description  of  the  spatially  devel¬ 
oping  twodimensional  mean  flow  be  recognized  and  Included  In  the  derivation  of 
the  vortlclty  conservation  equation  within  the  critical  layer  of  spatially 
amplifying  nonlinear  Instabilities.  Solutions  exhibiting  spatial  amplification, 
fixed  temporal  frequency,  but  gradually  varying  wavelength,  may  be  gleaned, 
albeit  laboriously,  for  such  an  equation  —  but  not  for  equation  (35),  which 
Is  predicated  on  the  assumption  of  streamwise  homogeneous  mean  flow.  Finally, 
in  connection  with  the  conclusion  Hi)  we  suggest  that  Robinson's  prediction 
of  a  fixed  value  for  the  growth  rate  ratios  of  the  harmonics  stems  from  the  local 
similarity  character  of  his  solution,  only  dependent  on  the  Instantaneous 
value  of  the  parameter  vij.  Experimentally  observed  growth  rate  ratios  scatter 
around  that  value  because  the  actual  harmonic  forcing  process,  described  in  a 
previous  paragraph,  does  not  obey  local  similarity. 

Aside  from  the  controversial  conclusions  discussed  above,  the  crucial  short¬ 
coming  of  Robinson's  model  resides  in  the  poor  comparison  between  Its  predic¬ 
tions  and  those  of  nonlinear  models  for  the  vorticity  distributions  associated 
with  nonlinear  instabilities  undergoing  saturation  and  subsequent  modulation. 

Due  to  the  unique  dependence  of  Robinson's  solution  upon  the  parameter 
the  comparison  becomes  especially  strained  at  those  post-saturation  stages  in 
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the  actual  modulatlonal  evolution  of  the  Instability  where  pj  has  the  same 
modulus,  but  opposite  sign  (e.g.  t«  19.15  vis  a  vis  t  ■  25.7  In  Figure  2); 
as  a  consequence,  the  modulation  itself  Is  not  predicted.  Whereas  that  pro¬ 
cess  Is  central  to  the  cyclic  mechanism  of  subharmonic  generation  and  vortex 
pairing  proposed  In  this  paper,  solutions  of  equation  (35),  subject  to  appro¬ 
priate  time  dependent  boundary  conditions,  are  prerequisite  to  a  realistic 
description  of  the  vortlclty  redistribution  within  the  critical  layer  of  a 
nonlinear  Instability.  Unfortunately,  such  solutions  are  necessarily  numerical 
and  rather  laborious,  as  they  require  the  parallel  consistent  determination 
of  the  matching,  slowly  varying,  outer  flow.  Unified  numerical  analysis  of 
the  critical  layer  and  the  outer  flow  by  a  fully  nonlinear  model  proves  more 
economical  in  practice  and,  therefore.  Is  adopted  here.  That  consideration 
notwithstanding,  the  inner/outer  layer  view  of  nonlinear  instabilities  is  use¬ 
ful  in  that  it  exhibits  the  vortlclty  redistribution  mechanism  operative  within 
the  critical  layer  as  well  as  the  limiting/dominant  influence  of  that  process 
upon  the  growth  of  the  instability  and  the  forcing  of  the  harmonics. 

The  conclusion,  repeatedly  suggested  above, that  only  nonlinear  models  —  which  do 
not  assume  either  local  similarity  or  self-preservation  of  the  wave  motion 
structure  —  realistically  describe  the  evolution  of  nonlinear  instabilities, 
is  also  supported  by  an  examination  of  the  approximate  model  suggested,  but  not 
applied,  by  Miura  and  Sato  (1978).  This  model  assumes  that  the  evolviny  struc¬ 
ture  of  a  temporal  instability  having  wavenumber  a  may  be  approximated  by 
(see  Equations  (3.11)  and  (3.13)  of  Miura  and  Sato  1978) 

e  ^(y.t)  =  2A(t)  f(y)  exp  [ie(t)  g(y) ]  (40) 
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and  that  the  attendant  harmonics  may  be  neglected.  Substitution  of  (40)  Into 
(29-a,b)  then  yields 

t 

U(y,t)  ■  U(y,0)  +  2a(fV)7  A*8dt 

o 

(see  equations  (3.20)  of  Mlura  and  Sato  1978)  as  well  as  two  nonlinear  ordinary 
differential  equations  for  A(t)  and  6(t),  which  follow  upon  formal  Integration 
of  the  real  and  Imaginary  parts  of  (29-b)  over  the  width  of  the  flow  (see 
equations  (3.14)  and  (3.17)  of  Mlura  and  Sato  1978).  The  stipulation  (40) 
clearly  implies  that  the  wave  motion  structure  retains  Invariant  spatial  depend¬ 
ences  f(y)  and  g(y)  for  the  amplitude  and  phase  functions*  respectively.  The 
results  of  Figure  5a  support  this  approximation  for  times  t  <  13.5,  prior  to 
saturation;  and.  In  fact,  an  Integration  of  the  equations  for  A(t)  and8(t)  — 
with  coefficients  based  on  functions  f(y)  and  g(y)  determined  by  the  eigen¬ 
mode  structure  at  t  *  0  —  yields  reasonable  predictions  for  the  presatura¬ 
tion  behavior  (Figure  2),  as  well  as  the  saturation  time  and  amplitude,  of 
the  test  case  instability  repeatedly  discussed  in  this  section.  Unfortu¬ 
nately, as  shown  in  the  next  section,  the  assumption  of  invariant  spatial  depend¬ 
ences,  collapses  at  later  times.  As  a  consequence,  our  solutions,  employing 
the  approximate  equations  derived  from  (40),  fail  to  reproduce  the  postsatura¬ 
tion.  oscillations  in  wave  amplitude  and  structure,  which  Miura  and  Sato 
claimed  to  extract,  but  did  not  calculate,  from  their  approximate  model.  Our 
central  interest  in  those  oscillations,  vividly  displcyed  by  Miura  and  Sato's 
own  solution  of  equations  (29-a,b)  as  well  as  by  more  elaborate  models  (e.g., 
ZabusRy  and  Deem  1971,  Patnaik  et  al  1976),  motivates  our  complete  reliance 
on  the  nonlinear  model  for  the  following  study  of  the  vortex  pairing  mechanism. 
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5.  NONLINEAR  MECHANISM  OF  SUBHARMONIC  EXCITATION  AND  AMPLIFICATION 

As  noted  in  previous  sections  of  the  paper,  measurements  of  time-averaged 
properties  as  well  as  fluctuation  spectra  in  forced  and  unforced,  transl- 
clonal  and  turbulent,  statistically  stationary,  spatially  developing  mixing 
layers  Indicate  that,  at  least  for  those  situations  where  the  forcing  fre¬ 
quency  Is  smaller  than  the  most  unstable  frequency  at  the  initial  station, 
the  specific  spectrum  of  disturbances  initially  coupled  to  the  flow  has  a 
vanishing  Influence  on  the  development  of  the  flow  downstream  of  the  station 
where  two  or  more  vortices  first  agglomerate  (Wygnanski  et  al  1979b,  Ho  and 
Huang  1980).  Following  that  event,  the  spatial  growth  rate  of  the  forced 
mixing  layer  thicxness  becomes  essentially  the  same  as  that  of  the  unforced 
layer  characterized  by  the  same  velocity  ratio  X*  (Uj  -  U2)/ (Uj  +  U2).  The 
finding  is  significant  In  that  it  associates  a  specific  physical  event  —  one 
ore  more  agglomerations  —  with  the  traditional  observation  that,  following 
some  initial  relaxation  distance,  unforced  mixing  layers  tend  to  grow  at  a  con¬ 
stant  rate,  which  depends  on  the  velocity  ratio  X,  but  Is  Insensitive  to 
the  specific  initial  conditions,  transition  tripping  device  and  environment  of 
small  amplitude  random  external  disturbances  encountered  in  any  one  experiment. 


Whereas  the  constant  growth  rate  of  the  self  preserving  unforced  mixing  layer  is 

paced  by  a  cyclic  process  of  vortex  pairing  over  statistically  averaged  lengths 

! 

which  scale  with  the  local  layer  thickness  6^  and  the  velocity  ratio  X,  the 
rapid  onset  of  the  same  dominant  process  in  forced  mixing  layers  strongly  sug¬ 
gest  that,  after  the  first  agglomeration,  the  history  of  the  large-scale 
structures/vortices  becomes  governed  by  a  self-reproducing  limit  cycle  of  non¬ 
linear  instabilities,  viz  the  cycle  already  described  in  Section  II  in  the 


T-i' 
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course  of  Introducing  and  justifying  the  nonlinear  model  used  In  the  present 
paper.  The  key  features  of  the  cycle »  which  we  seek  to  analyze  and  demonstrate 
here*  are:  1)  the  self-induced  modulation  of  the  nonlinear  Instability 
responsible  for  an  agglomeration  of  vortices  In  the  mixing  layer  combines 
with  the  concurrent  distortion  of  the  mean  flow  to  provide  Intrinsic  excita¬ 
tion  of  subharmonic,  spatially  amplifying,.  Instabilities  having  about  half 
the  frequency  of,  but  nearly  the  same  phase  velocity  as, the  self-modulating 
wave,  hereafter  identified  as  the  fundamental;  11)  the  Initial  finite  ampli¬ 
tude, phase  and  structure  of  the  subharmonlc  wave  motion  produced  by  the 
modulation  of  the  fundamental- Independent  of  either  local  or  previous  random 
external  forcing  at  the  appropriate  frequency  —  are  conducive  to  nonlinear 
amplification  of  the  subharmonic,  which  draws  energy  from  the  fundamental  as 
well  as  from  the  mean  flow;  ill)  because  of  the  specific  phase  relationship 
between  the  fundamental  and  the  amplifying  subharmonic,  pairs  of  vortices, pre¬ 
viously  generated  by  the  fundamental,  merge  with  the  consequence  that  the  funda¬ 
mental  becomes  essentially  extinct,  while  the  subharmonlc  becomes  dominant; 

1v)  due  to  this  predominance  and  to  the  self-limiting  nature  of  the  subharmonlc 
growth  —  necessarily  bounded  by  the  finite  amount  of  vortlcity  available  to 
the  wave  within  its  fixed  wavelength  —  specific  structures  of  the  mean  flow 
and  of  the  wave  motion  evolve  at  subharmonic  saturation,  viz  structures  which' 
are  conducive  to  the  self-preserving  repetition  of  the  instability  cycle  because 
they  reproduce,  except  for  a  doubling  of  the  length  scales,  the  state  of 
motion  prevailing  at  the  earlier  time  of  fundamental  saturation. 

Representative  time  dependent  structures  of  the  mean  flow  and  of  the  wave  motion 
which  evolve  upon  saturation  of  a  dominant  finite  amplitude  instability,  are 
illustrated  in  Figures  2,  6,  8  and  9,  which  display  the  long  time  evolution 
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(t  >  13.5)  of  the  specific,  temporally  developing, shear  layer  discussed  in 
Section  4,  viz  the  flow,  bounded  by  rigid  plates  at  y  =  ±  3.875,  which 
evolves  from  the  initial  mean  velocity  field  U(y,Q)  =  tanh2y  and  the  super¬ 
posed  unstable  eigenmode  with  a  =  0.8892,  p(0)  =  0  +  -.37841,  e  =0.1, 

<j>l(0,0)  =0.1,  under  the  assumption  that  the  motion  remains  at  all  times 
spatially  periodic  with  wavelength  (2II/a).  Clearly  the  vorticity  distribu¬ 
tion  achieved  at  the  time  when  the  instability  first  saturates  is  not  in 
dynamic  equilibrium  for  the  mean  as  well  as  the  wave  motions  undergo  subse¬ 
quent  structure/ amplitude  modulations  with  period  T  «  17,  about  2.4  times 
as  long  as  the  period  Tfa7.06  of  the-  considered  instability  propagating 
with  phase  velocity  Cp  =  1  in  the  mean  flow  of  initial  velocity  profile 
U(y,0)  =  (1  +  tanh2y).  Although  periodic,  the  noted  modulations  are  not 

strictly  sinusoidal:  for  example.  Figure  2  shows  that  the  stationary  value 

/ 

of  the  y-derivative  of  the  phase  function  on  the  axis  y  *  0/ 

0,y(O,t)  =  { (3/3y)  tan"1(Im[<|,1(y,t)]/Re[<|>1(y,t)3)}ysp 

occurring  at  t  s  19.15  precedes  the  zero  of  the  same  function  (and  the  concur¬ 
rent  minimum  of  the  disturbance  energy)  occurring  at  t=22  by  a  time  interval 
At  «  2.85,  /which  corresponds  to  a  period  ^  «  11.4.  Thus,  the  modulation 
encompasses  a  band  of  frequencies  centered  around  4n  /(Tmu  +  T^)  =  2n/Tma, 
where  Tm,s  14.2.  Since  the  period  T  =  2  Tf  almost  exactly,  the  frequency 
spectrum  of  the  self-modulating  wave  motion  in  the  flow  characterized  by  the 
initial  mean  velocity  'Ufy^JV®  (1  +  tanh2y)  would  then  exhibit  -  in  additon 
to  discrete  peaks  at  the  frequencies  (2nII/Tf)  of  the  dominant  wave  (n  =  1) 
and  its  hairmonics  (n>  1)  -  a  band  of  width  ±  0.2  (n/T^.)  centered  at  the 
frequency  (n/Tf)  of  the  1/2  subharmonic.  Interestingly  enough,  the  excitation 
of  just  such  a  frequency  band  is  displayed  by  the  spectral  measurements  made 
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downstream  of  the  station  where  the  fundamental  saturates  in  a  transitional  mixing 
layer  forced  at  the  most  unstable  frequency  (see  Mlksad  1972,  Figure  8,  x>  11.5) 

Since  an  approximate  doubling  of  the  wavelength  relative  to  that  of  the  funda¬ 
mental  is  Inherent  to  the  subharmonics  observed  In  the  experiments,  but  not  to 
the  modulation  of  the  time -dev el  oping  shear  layer  considered  In  the  calculations, 
some  reflection  on  the  stipulated  equivalence  between  temporally  and  spatially 
developing  flows  is  in  order. 

For  simplicity  of  discussion,  the  relatively  narrow  band  of  modulation  frequencies 
exhibited  in  Figures  2,  6,  8  and  9  is  compressed  into  the  central  1/2  subharmonic 
with  strictly  trigonometric  time  dependence.  Accordingly,  the  structure  of 
the  temporally  developing  flow  at  times  £  *  (t  -  ts)  >  0  (where  t$  =  13.5  corre¬ 
sponds  to  saturation  of  the  fundamental)  is  approximated  by 

U(y,0  *  U(o)(y)  +  e[U(l)(y)  *  c.c.)  (41-a) 


eax  =  e^1(0',t)[4i1(y,t)/<til(0:,t)]eiaX  = 
e*l<°%1(y,t)  M  .  *{*/«>  .  Wl(1>  +  c.c.Dx  l«-b) 

{?J0)(y)  +  i*,<0)(y)  (y)  «1l,t/Tf  +  c.c.]  ♦  Il*1<1)(y)«1,,t/Tf  ^ 


lax 


the  complex  functions/quantities  U^(y),  (y),  (y),  m  =  0,1 

being  defined  by 


,(m)  m  2  T 


(y)  =  (2em  Tf)_1  /  f  U(y,t)  ^imiTt/Tf  dt 


(42-a) 
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(m) 


*  (2Tf )_1  ff  4,^0, t)  «,mIr'/Tf  dt 
0 

=  (2Tf)‘X  /Tf  dt 


(42-b) 


(42 -c) 


S  (m)(y)  -  (2Tf )-1  /Tf  1m  tfyy.t)]  eimt/Tf  dt 


(42-d) 


Calculations,  as  well  as  inspection  of  Figures  2,  6,  8  and  9 indicate  that 
|U(1>|  =  o(l),|d1(1)|<|  d^l,  I?ic(1)h  oU|¥1<°)|  and|?1'1)|<|?1>)|for 
the  problem  on  hand.  Thus,  if  only  terms  of  first  order  are  retained  in  the 
modulation  function/quantities  characterized  by  the  superscript  m  -  1,  the 
stream  function  $  ^  descriptive  of  the  time  developing  flow  in  the  domain 
T  £  0  may  be  cast  i n  the  form 

${t)(x,y,t)  -  Y(y)  +  £<i>(x,y,t)  =  /  U(o)  (y)  dy  +  |  {2  /u^i(y)dyeiBrt/2  + 

yw  yw  ‘  _ 


E  f,W)  (y)ei[ax-{1-j)8rI]+c.c.) 

2j=i,2,3  1 


where  6r  =  2-nrTf “  and 

f<1/2)(y)  -  4,'11*  *lc<°>(y)  *  .  *ls(1)*(y) 


(44-a) 


fj'11  (y)  ■  *1<0>  4ic(0)(y) 


(44-b) 


f!(3/2)(y)  -  *i(1)  *ic(0)(y)  + 1  *i(0)  ?is(1)  (y) 


(44-c) 
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This  approximate  representation,  together  with  the  results  presented  In 
Section  4,  provides  the  context  for  examining  the  structure  of  the 
laboratory  mixing  layer  and  the  applicability  of  the  equivalence  in  the 

i 


spatial  domains  x^  <  x^s  and  x^  >  x^.which,  respectively,  correspond  to 
t  <  o  and  "t  >  o. 


The  equivalence  of  temporally  and  spatially  developing  flows  presumes  that  the 
approximate  relationship  between  linear  temporal  and  spatial  instabilities  estab¬ 
lished  by  Gaster  (1962)  applies  as  well  to  the  nonlinear  flows  considered  here. 
Accordingly,  the  stream  function (x,y,t),  which  describes  the  evolving 
structure  of  the  theoretical,  monochromatically  forced,  unsteady  shear  layer 
viewed  in  coordinates  (x,y)  moving  at  the  phase/group  velocity  c  sc  =(6r/a)=l 
of  the  forcing  instability  with  real  wavenumber  a  and  complex  frequency  (3f+  ipj) 
and  the  stream  function  ^^(x^,  y^,  t^),  which  describes  the  structure  of  the 
experimental,  statistically  stationary,  mixing  layer  viewed  in  laboratory  fixed 
coordinates  (x^.y^)  and  subject  to  the  spatially  amplifying  instability  with 


complex  wavenumber  a^+io^  oil  +  .i(3j/Br)]  real  frequency  8 and  nearly  coinci¬ 
dent  phase/group  velocities  cps cg*  8/a^l,  are  related  by 

$^(rix,  y.B^)  ■  (ar  x  -  Bt. ,  yA,  or  x^)  -  /*!  (p/ar)  dy£  (45) 


when  the  zero  of  t.  is  chosen  so  that 
i 


■I* 


ax,  y,  0)  =  ^S^(-Bt  ,  y  ,  0)  -  /(3/a  )dy 


l  l  '  '  —  r-  ’ i 

~yAw 

However,  (45)  has  a  limited  domain  of  validity,  as  discussed  below. 

A  first  limitation  to  the  equivalence  becomes  manifest  upon  inspection  of  the 
eigenfunctions  and  eigenvalues  calculated  in  the  linear  limit:  for  the  initial 
mean  velocity  profile  U  =  (1  +  tanh2y)  considered  here,  only  spatial  and  tem¬ 
poral  instabilities  with  dimensionless  frequency  Re(B)  >0.8  approximately 
satisfy  the  conditions  on  wavenumbers,  wave  velocities  and  wave  motion  struc¬ 


3 

i 

'1 

1 

1 

j 


tures  prerequisite  to  (45).  Spatial  instabilities  which  comply  with  the 
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equivalence  conditions  at  small  amplitudes  (xfc  *  o)  continue  to  do  so  In  the  | 

process  of  nonlinear  growth  and  saturation  at  stations  <  x^s.  In  that  domain  5 
the  phase  and  group  velocities  of  the  spatially  amplifying  waves  remain  essen-  | 

:i 

tlally  invariant  because  the  effective  wavenumber  ae(xA)  =  [<*r  Uy(0.0)/Uy(xif0)I  j 
increases  with  xfl.  The  velocity  field  within  the  critical  layer  may  then  be  J 

approximated  by  | 

I  .4 

I  I 

ufc  =  (B/ar)  +  Uy(x*,o)y*  «  1  +  V  V  0)y* 

v^  *  a(xA)  or  s1n(ar  x&-  Bt^)  +  o(a*/*)  | 

and  the  vorticlty  field  by  ] 

v!f  <*><,,.  y,,  tt)  ■  u  (0.0)  /  +  afXj)  a,*5*  (xr  Y.  n)  . ] 

(s)  i 

where  n  =  (af  x^  -fJt^)  and  Qj  '  denotes  a  function  assumed  of  order  unity  1 

toq ether  with  its  derivative  with  respect  to  the  stretched  coordinate 

i 

j 

Y  =  lV  V  °)/a(xA)l%  y*,  =  (x£)l  } 


across  the  critical  layer.  To  a  first  approximation  the  function  8.  ^  •  i 

which  describes  the  concurrent  nonlinear  spatial  evolutions  of  the  (time 
averaged)  mean  flow  and  the  wave  motion  within  the  critical  layer  —  satisfies 

I  H 

the  equation 


+  Y8.<s)  +sinnniJs)  +  u1Ini(s)-(.1+ w2/u1)-|-niY(s)]  =  o 

2 

in  boundary  layer  coordinates  (X,  Y,  r\)  defined  by  dX  =  (ar  a/36)dxA, 
Y  =  (y£/6),  n  =(ar  x^  -  Btjjh  the  quantities 


y^X)  =  Id  loga(x)/dX],  u2  “  "  ld  log  U^(x»0)/dX] 


L. 
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being  of  comparable  magnitude  <1  In  the  domain  of  nonlinear  growth.  Equation  j 
(46)  Is  formally  Identical  to  the  equation  (35)  previously  derived  for  tempor-  j 

ally  developing  shear  layers,  and  so  are  the  attendant  Initial  conditions  that  J 

84^(0,  Y, 11)  be  periodic  In  n  as  well  as  the  boundary  conditions  that  for  J 

Y  ■+*••,  the  velocities  and  vortlclty  associated  with  0^  S.(X»  Y »yi)  match  3 

I 

those  associated  with  a  mean  flow  slowly  varying  In  x  P^s  *  superposed  spa-  j 
tlally  amplifying  weakly  nonlinear  wave  having  fixed  frequency  8  and  wavenumber  j 
oA  that  varies  slowly  with  x-  As  a  result  a/s)(X,  Y,n)  -  ^(T.Y.E)  for  xA<  xAs,| 
t  .<  0,  and  equivalence  Is  preserved  In  those  domains.  •  j 

Ihe  equivalence  becomes  suspect  for  xA>  xAs  because  the  structure  of  the  j 

laboratory  mixing  layer  predicted  by  (43)  and  (45)  • 


t>1 


ID  ■  (y.)  *  * 


u^(y,)  dy '  W»)  +' 


W) 


f’1(1)(yj,)*t(or  **  -«*)  (ye)e,(6/2)(te  *  *e  '  e V  +  c.c.i 

loses  strict  time  periodicity  at  the  forcing  frequency  8,  while  the  theoretical 

*  4  \ 

time  developing  shear  layer  retains  spatial  periodicity  with  the  Initial 

» 

wavenumber  o.  In  addition  to  the  fundamental  —  with  fixed  amplitude,  struc¬ 
ture  f/^(yA)  ar,d  'requency  p  —  (47)  Indicates  the  presence  of  the  1/2  and 
3/2  harmonics,  with  structures  subject  to  a  slow  time  modulation  at  the 
frequency  8/2.  Linear  analysis  readily  shows  that,  for  the  prevailing  mean 
velocity  profile  [(p/cy)  +  (y^) ]  »  the  1/2  harmonic  constitutes  an  admissi¬ 

ble  temporal  as  well  as  spatial  instability  with  nearly  equal  phase  and  group 
velocities  cp=jcg  s(p/ar)  practically  coincident  with  those  of  the 
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dominant  fundamental,  but  the  3/2  harmonic  does  not.  The  possibility  then 
arises  that  a  time-modulated,  or  Intermittent,  excitation  of  the  1/2  subhar- 

=j 

monlc  Instability  occur  at  appropriate  downstream  stations  in  a  spatially 
developing  mixing  layer  dominated  by  a  finite  amplitude  large  scale  Instability 
as  the  counterpart  to  the  temporal  modulation  of  the  fundamental  wave  motion.  'j 

j 

structure  In  the  equivalent  temporally  developing  shear  layer.  Should  the  | 

subharmonic  Instability  become  amplified,  the  equivalence  between  the  two  \ 

flows  would  clearly  be  broken.  j 

Linear  analysis  falls  to  Indicate  whether  and  when  the  1/2  subharmonic  Is  ampli-  | 

fled  because  the  stream  function  of  the  associated,  slowly  modulated  wave 

motion  with  real  part  j 

. .  ■  "i 

IRe(f^1/2V2fy  dyAlcos  §  1*1  "  *£5)  +  Im(f^1/2))sin  f  (t^-  t£$)  (48-a)  j 

and  Imaginary  part 

lm(f1(1/2))  cos  |  l,t,  -  tls)  -  tRe(f1(1/2))  -  2 f*  U(1)  (48-b) 

carries  little  resemblance  to  that  of  a  linear  wave  at  any  time  0  <  (t^  -  t^  ) 

•j 

<  4tt/0.  This  structural  difference  also  prevents  the  adoption  of  either  Kelly's 

3 

(196?)  theory  or  its  extensions  in  terms  of  weakly  nonlinear  multiple  scale 
models  of  secondary  instabilities  (e.g.,  Nayfeh  and  Bozatli  1979)  to  determine 

i 

j 

the  paranrteric  amplifier  action,  if  any,  of  the  coexistent  nondisperslve  finite  j 
amplitude  fundamental.  Seemingly,  the  only  resort  then  resides  in  nonlinear  i 

numerical  calculations,  which  in  view  of  the  nondisperslve  nature  of  the  funds-  | 
mental  and  the  subharmonic  at  the  considered  stage  of  flow  development,  may  be 
based  on  the  heuristical ly  equivalent,  time-developing,  nonlinear  flow  model  of  1 


equations  (29,a,b,c,). 
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Numerical  experiments.  Involving  the  Integration  of  (29-a,b,c)  subject  to 
the  distinct  sets  of  Initial  conditions  which  accrue  when  the  mean  velocity 
field  and  the  fundamental  wave  motion  displayed  In  (47)  are  combined  with  the 
subharmonic  wave  motions  indicated  by  (48-a,b)  for  distinct  tlnies  t^^  Indicate 
that,  subject  to  u  variable  time  delay  generally  shorter  than  (2/p),  nonlinear 
effects  cause  the  subharmonlc  to  undergo  a  strong  nearly  exponential  amplifi¬ 
cation  with  characteristic  time  comparable  to  (2/0)  for  those  Initial  conditions 
which  pertain  to  times  t^  In  the  interval  roughly  bounded  by  - (tt/3)  <  (t^  -  t^s) 

<  (it/ p).*  Since  those  Initial  conditions  recur  with  a  time  period  (4ir/0) 

according  to  (48-a,b),  strong,  but  intermittent,  subhanmonic  growth  is  indicated 
for  the  region  of  the  laboratory  mixing  layer  which  extends  approximately  (2^/0^,) 
downstream  of  the  station  =  x£$  where  the  forcing  fundamental  saturates.  Such 
a  region  of  intermittent  nonlinear  subharmonic  excitation/amplification  --  Inher¬ 
ent  to  the  evolution  of  the  fundamental  and  Independent  of  external  forcing  at 
the  appropriate  frequency  --  has  been  observed  by  Miksad  (1972,  p.  714)  in 
transitional  mixing  layers.  Its  occurrence  and  interpretation  are  quite  signifi¬ 
cant  in  our  view  of  mixing  layer  dynamics  because  they  provide  the  first  two  ele¬ 
ments  of  the  self-reproducing  nonlinear  cycle  of  large  scale  instabilities  which 
we  seek  to  demonstrate. 

e  nonlinear  nature  of  the  subharmonic  excitation/amplification  process  requires 
that  it  be  analyzed  in  the  context  of  the  concurrent  changes  in  mean  flow 

and  fundamental  wave  motion.  The  intermittent  activity  over  the  noted  distance 
t  ?tream  of  x„s  admittedly  pre/ents  the  exact  simulation  of  that  region  of  the 

*A  reexamination  of  the  evolving  structure  of  the  fundamental  wave  motion  prior 
to  saturation  (see  Section  4)  shows  that^the  stream  function  is  adequately 
described  by  the  approximation  (43)  for  T  >-(ir/0). 
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laboratory  flow  by  a  single,  equivalent,  time-developing  shear  layer  having 
the  streamwlse  periodic  structure  attendant  to  superposed  sinusoidal  funda¬ 
mental  and  subharmonic  waves.  However,  intuition  as  well  as  data  (e.g. ,  Figure 
3)  suggest  that  only  a  modest  amount  of  energy  accrues  to  the  subharmonic  in 
the  region  of  intermittent  activity.  An  approximate  simulation  then  seems 
plausible,  especially  if  attention  is  focused  on  the  flow  structure  which 
ultimately  evolves  upon  finite  amplitude  growth  and  saturation  of  the  subhar¬ 
monic.  In  that  connection  a  station  xAo  Is  chosen  within  the  region  in  question 
downstream  of  xa$;  non-intermlttent  activity  is  assumed  to  begin  at  that  sta¬ 
tion,  and  effective  initial  structures  are  prescribed  there  for  the  time- 
averaged  velocity  profile  as  well  as  for  the  fundamental  and  the  subharmonic 
wave  motions.  Given  the  small  energy  accrued  to  the  subharmonic,  the  effective 
initial  structures  may  not  differ  substantively,  from  those  indicated  by  (47); 
accordingly,  a  parametric  study  need  only  consider  moderate  variations  in  the 
partition  of  energy  among  the  modes  and  in  the  structure  of  the  subharmonic  wave 
motion  from  its  time-average  during  the  periods  of  activity.  Limitpd  numerical 
experiments  indicate  that  such  variations  have  little  influence  on  the  ultimate 
development  of  the  flow.  In  the  absence  of  a  cogent  criterion  to  guide  the 
choice,  only  the  representative  results  of  a  single  approximate  simulation  are 
therefore  reportedand  discussed  below;  the  pertaining  initial  conditions  are: 
i)  xAo  =  x  +  (Tr/af);  (ii)  structures  of  the  mean  velocity  profile  and  of 
the  fundamental  wave  motion  as  displayed  in  (47);  iii)  structure  of  the  sub¬ 
harmonic  motion  as  indicated  by  (4/)  for  (t^  -  t^  )  a  0 

I 

The  structural  evolution  of  the  equivalent  temporally  developing  if low,  which 
ensues  upon  excitation  of  the  subharmonic,  is  illustrated  in  Figures  10,  11  and 
12.  Inspection  of  these  figures  elicits  several  comments.  Clearly  the  subhar- 
monic  amplification  and  saturation  constitute  strongly  nonlinear  processes 


-44- 


NYU/DAS  81-06 


accompanied  by  large  changes  In  the  structures  of  the  mean  velocity  profile 
as  well  as  the  fundamental  and  the  subharmonlc  wave  motions.  The  salient  conse¬ 
quences  Include:  1)  the  binary  merging  of  the  vortlcity  concentrations  generated 
by  the  preexisting  fundamental  which,  in  fact,  transfers  most  of  its  energy  to 

5 

the  subharmonic  as  the  latter  approaches  saturation  (Figure  12);  11)  an  increase  I 
In  the  mean  flow  vortlcity  thickness  6^  which,  at  subharmonic  saturation,  becomes 

very  nearly  double  that  prevailing  at  fundamental  saturation  (Figures  8  and  10);  } 

-J 

ill)  an  evolution  in  the  subharmonic  wave  motion  which,  consistent  with  flow  self-  j 

preservation,  acquires  at  saturation  a  structure  and  an  energy  which  are,  respec-  j 

1  } 

tively,  equal  and  double  their  counterparts  in  the  fundamental  wave  motion  at  its  -i 

I 

saturation  (Figures  9,11  and  12);  iv)  a  self-induced  modulation  of  the  mean  and 

subharmonic  wave  motions  which  follows  subh-rmonic  saturation  and,  consistent  j 

J 

with  the  flow  structure  and  the  time/length  scales  attendant  to  that  event,  causes  i 
first  the  intrinsic  excitation  of  a  1/4  subharmonic  wave,  and  then  its  evolution  l 

3 

with  a  life  cycle  dynamically  identical  to  that  of  the  subharmonic.  In  combi-  | 
nation  these  processes  complete  the  demonstration  of  the  proposition  that  a 
single  self -reproducing  cyclic  process  of  1/2  subharmonic  excitation,  amplifi- 

■i 

cation  and  quenching  controls  the  evolution  of  large  scale  coherent  structures 

j 

as  well  as  mean  flow  in  spatially  developing  mixing  layers.  1 


Disagreement  between  relevant  experimental  data  prevents  a  definitive  quantitative 
evaluation  of  the  analytical  results  for  the  nonlinear  wave  cycle.  The  sub¬ 
harmonic  amplification  rate  displayed  in  Figure  12  'is  measurably  larger  than 
that  reported  by  Miksad  (1972)  and  reproduced  in  Figure  3,  but  agrees  nicely 
with  that  observed  by  Ho  and  Huang  (1980,  see  Figure  17)  in  mixing  layers 
forced  at  frequencies  slightly  below  that  of  maximum  amplification.  The  apparent 
contradiction  may  be  facility  related.  Both  experiments  were  carried  out  in 
channels  with  relatively  small  cross-sections,  viz  12.5  x  13.5  cm  for  an 
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estimated  vorticlty  thickness  6^  «  3  cm  at  subharmonic  saturation  In  Miksad's 
case  and  10  x  10  cm  for  6w  «  1.*  cm  In  Ho  and  Huang's  case.  As  Mlksad  points 
out,  cross-stream  boundaries  then  exert  spanwlse  constraints  on  the  three 
dimensional  secondary  instabilities  with  large  spanwlse  wavelength  (see  Browand 
and  Troutt  1980,  Figure  6)  which  arise  In  the  course  of  subharmonic  amplifica¬ 
tion,  and  cause  skewing  as  well  as  apparent  branching  of  the  large  scale  struc¬ 
tures  (see  Miksad  1972  Figure  24  as  well  as  Browand  and  Troutt  1980).  Accord¬ 
ing  to  Miksad  "the  onset  of  spanwlse  activity  coincides  with  a  noticeable 

decay  of  fundamental  mode  energy  along  the  center- plane";  however,"  u'  _  _ 

■  •in*  s  * 

wave- front  measurements  Indicate  that  this  may  be  due  to  a  spanwise  redistribu¬ 
tion  of  energy  and  not  to  a  transfer  of  energy  to  other  scales  of  motion."  Thus, 
It  is  possible  that  the  subharmonic  amplification  rate  measured  by  Miksad  in 
the  center-plane  of  his  facility,  with  a  ratio  (d/6u)  “  4  between  channel  width 
d  and  mixing  layer  thickness  6^,  may  not  represent  either  the  spanwise  average 
value  or  the  value  appropriate  to  a  flow  free  of  lateral  constraints \  by 
contrast,  the  rate  observed  by  Ho  and  Huang  with  (d/6^)  =  8  may  not  be  subject  to 
such  adverse  limitations.  Admittedly  this  is  only  a  conjecture.  However,  a 
stronger  argument  ma;/  be  made  in  support  of  the  calculated  results  and  Ho  and 
Huang's  measurements  by  examing  the  average  spatial  growth  rate  of  the  mixing 
layer  over  the  interval  between  fundamental  and  subharmonic  saturation.  For 
the  mean  velocity  profiles  shown  in  Figure  10  and  the  considered  velocity  ratio 
X  -  1  we  obtain 

(dQ/  dx)  =  (de/dt)  » (l/4)(d6w  /dt)  *  .03 

a  result  which  scales  in  the  lower  range  of  the  experimental  observations  reported 
by  Brown  and  Roshko  (1974)  and  Browand  and  Latigo  (1979),  if  the  growth  rate  of 
the  integral  (momentum)  thickness  0  is  assumed  to  depend  linearly  on  X.  Thus, the 
the  most  commonly  measured  gross  feature  of  mixing  layer  development  is  fairly 

reproduced  by  the  model  and  calculation  presented  here. 
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6.  DISCUSSION  | 

The  results  of  section  5  Indicate  that  the  growth  of  spatially  developing 
turbulent  mixing  layers  Is  controlled  by  a  self-reproducing  deterministic,  non¬ 
linear  wave  mechanism  whereby  sequentially  subharmonic,  large-scale,  finite 
amplitude  Instabilities  are  caused  to  follow  a  well  defined  excitation/ ampli¬ 
fication/quenching  life  cycle,  which  renders  them  sequentially  dominant.  The 
key  to  the  mechanism  resides  in  the  strongly  nonlinear  nature  of  the  Inter¬ 
actions  (between,  and  the  concurrent  structural  evolutions  of,  the  mean  flow 
and  the  locally  dominant  wave  motions.  The  salient  consequences  of  those 
nonlinearities  Include:  i)  the  achievement  of  negligible  dispersion  between 
a  dominant  instability  at/near  saturation  and  Its  1/2  subharmonic,  provided 
the  frequency  of  the  fundamental  falls  In  a  range  to  be  discussed  below; 

11)  the  onset  of  non-stat Ionary  mean  and  wave  motions  upon  saturation  of  the 
dominant  fundamental,  and  the  consequent  intermittent  excitation  of  1/2  and 
3/2  harmonic  wav '  motions  having  well  defined  structures;  iii)  the  approximate 
doubling  of  the  mean  flow  thickness  and  the  near,  complete  draining  of  the 
fundamental  energy,  which  accompany  the  amplification/saturation  of  the  unstable 
subharmonic  ^and  permit  the  self- preserving  reproduction  of  the  (subharmonic) 
instability  excitation/amplification/quenching  processes  in  sequential  cycles 
distinguished  only  by  a  sequential  doubling  of  the  length  of  scales  involved. 

The  deterministic  and  skeletal  nature  of  the  proposed  cyclical  instability  mech¬ 
anism  --  ascribing  a  significant  role  to  a  discrete,  two  frequency,  spectrum 
of  two  dimensional  waves  in  each  cycle  —  has  distinct  implications  with  regard 
to  the  dynamics  of  turbulent  mixing  layers  and  their  amenability  to  control. 

The  quantitative  pursuit  of  these  implications  is  beyond  the  scope  of  the  pre¬ 
sent  paper;  however .some  relevant  qualitative  trends  may  be  gleaned  from  the 
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available  results.  The  first  comment  In  this  connection  concerns  the  expected 
behavior  of  transitional  mixing  layers  monochromatlcally  forced  at  unstable 
frequencies  8  other  than  the  ^  of  maximum  amplification  considered  in  our 
calculations.  Since  the  growth  of  any  instability  is  limited  by  the  finite 
amount  of  vorticity  available  to  it  within  a  wavelength,  one  anticipates  that. 

In  accord  with  the  observations  of  Miksady (1972,  p.  712),  the  saturation  ampli¬ 
tude  of  a  forced  disturbance  Increase  as  its  frequency  decreases.  The  mean 
flow  distortion  (Increase  In  6^),  which  accompanies  the  growth  of  the  Instability 
and  is  geared  to  its  saturation  amplitude,  is  also  expected  to  increase  as 
the  frequency  decreases.  As  a  result,  conditions  conducive  to  the  excitation 
of  a  1/2  subharmonic,  having  nearly  the  same  propagation  velocity  as  the 
forcing  fundamental,  and  the  attendant  initiation  of  a  nonlinear  subharmonic 
generation  cycle.may  be  anticipated  for  a  finite  range  of /forcing  frequencies 
about  6m.  The  experiments  of  Ho  and  Huang  confirm  this  expectation  for 
(Bm/2)  <8  <  Bm*»  the  range  of  admissible  frequency  8  >  Bm  remains  to  be  deter¬ 
mined  by  either  calculation  or  experiment. 

For  a  specified  initial  mean  flow,  forcing  at  different  frequencies  in  the 
aforenoted  range  results  in  initial  mixing  layer  responses,  hereafter  identified 
as  mode  1,  which  are  qualitatively  identical,  but  quantitatively  different 
because  each  frequency  possesses  a  specific  admittance,  amplification  rate, 
saturation  amplitude  and  rate  of  energy  transfer  to  its  subharmonic.  However, 
the  frequency  sensitivity  should  largely  vanish  upon  the  first  vortex  pairing 
because  the  subharmonic  responsible  for  that  event,  saturated  at  a  self-limiting 
amplitude  consistent  with  its  wavelength,  combines  with  the  attendant  commen¬ 
surate  distortion  of  the  mean  velocity  profile  to  yield  a  flow  structure  which, 
in  coordinates  scaled  to  the  local  mixing  layer  thickness,  is  independent  of  the 
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Initial  conditions.  If  the  nonlinear  Instability  cycle  proposed  here  Is 
operative,  the  development  of  the  flow,  downstream  of  the  station  where  vortex 
pairing  first  occurs,  should  then  become  self-preserving  and  Independent  of 
the  specific  transition  process,  whether  forced  or  natural.  This  .self-pre- 
serving  behavior,  which  is  the  cornerstone  of  the  proposed  model.  Is  manifested 
by  the  measurements  of  Ho  and  Huang  (1980)  In  forced  flows,  as  well  as  by 
Broward  and  Latigo's  (1979)  observation  that.  In  naturally  developing  flows 
subject  to  distinct  Initial  conditions,  the  mean  velocity  distributions  achieve 
an  equilibrium  shape  In  coordinates  scaled  to  the  local  mixing  layer  thickness 
after  two  vortex  pairings,  and  the  distributions  of  rms  longitudinal  fluctuations 
as  well  as  Reynolds  stress  do  so  prior  to  the  third  pairing.  Consistency  with 
the  proposed  mechanism  Is  ihus  indicated,  but  not  conclusive  support.  Clearly  , 
the  suggested  predominance  of  a  deterministic  two  dimensional  instability 
cycle  involving  a  discrete  two  frequency  spectrum  must  be  reconciled  with  the 
experimental  evidence  that  mixing  layer  turbulence  is  random,  three  dimensional 
and  broad  spectrum.  This  reconciliation  is  pursued  below  following  some  addi¬ 
tional  discussion  of  initial  mixing  layer  responses. 

Ho  and  Huang  (1980)  point  out  that  the  initial  response  changes  at  low  forcing 
frequencies  £  <  (8m/2).*  Different  modes,  designated  by  n  «  2,3,4,  arise  for 
different  frequency  ranges  ^/(n+1)  <3<8m/n.  Due  to  the  receptivity  and  ampli¬ 
fication  characteristics  of  the  flow.the  nth  harmonic  of  the  forcing  frequency 
dominates  the  initial  development  in  mode  n.  In  the  higher  modes  (i.e.  n  =3,4) 
the  (n-m)>l  harmonics  also  appear  and  initially  dominate  over  the  forcing 

Attention  is  here  restricted  to  small  amplitude  forcing,  which  has  a  small 
influence  on  the  spatial  evolution  of  the  instabilities  compared  to  their 
interactions  with  the  mean  flow  and  with  each  other.  The  collective  inter¬ 
actions,  which  sometime  arise  under  large  amplitude  forcing  (Ho  and  Nosseir 

1981),  require  separate  analysis. 
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frequency;  however,  their  energies  and  amplification  rates  remain  distinctly 
smaller  than  those  accruing  to  the  n**1  harmonic.  This  "response"  frequency 
then  amplifies  and  saturates  along  the  lines  discussed  in  Section  4,  with 
little  Interference  by  the  other  instabilities  simultaneously  present  in 
the  flow;  as  it  so  evolves,  it  generates  vorticity  concentrations  as  well 
as  mean  flow  distortions,  which  tend  to  eliminate  the  initial  differences  be-; 
tween  the  phase  velocities  of  the  various  instabilities. 

According  to  the  considerations  set  forth  in  a  previous  paragraph,  the  flow 
distortion  attendant  to  the  saturation  of  the  response  frequency  eliminates 
the  dispersivity  of  that  wave  relative  to  its  1/2  subharmonic,  which  coincides 
with  the  forcing  frequency  in  mode  2.  Binary  vortex  merging  then  ensues  in 
mode  2,  much  along  the  lines  repeatedly  discussed  in  this  paper,  but  without 
the  delay  required  for  intrinsic  excitation  of  the  subharmonic,  which  is  now 
forced  and  amplified  ab  initio.  the  forced  subharmonic  saturates  at  its  self- 
limiting  amplitude,  conditions  conducive  to  self-preserving  development  of  the 
flow  are  recovered  upon  the  first  vortex  pairing.  The  ensuing  development 
must  be  the  same  as  that  of  mixing  layers  initially  responding  in  mode  1,  and 
that  is  the  case  according  to  the  observations  of  Ho  and  Huang  (1980).  ( 

In  mode  3  the  saturation  of  the  response  frequency  is  also  followed  by  the 
continued  growth  of  its  1/2  su'jharmonic,  forced  and  amplified  ab  initio  in 
the  mixing  layer.  However,  that  process  does  not  proceed  to  self-limiting 
saturation  as  in  mode  2,  because  the  attendant  mean  flow  distortion  eliminates 
the  dispersivity  of  the  forcing  frequency  relative  to  the  response  frequency 
and  the  1/2  subnarmonic.  A  three-wave  nonlinear  interaction  then  arises, 
which  results  in  the  merging  of  every  triad  of  vortices  generated  by  the 
response  frequency.  Since  the  rate  of  energy  transfer  among  the  three  inter¬ 
acting^  waves! depends  upon  their  phase  relationship,  the  pattern  of  merging 
varies  from  event  to  event:  in  most  cases  two  members  of  the  vortex  triad 
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merge  first  Into  a  new  vortex,  which  then  coalesces  with  the  third  member; 
occasionally,  all  three  members | merge  at  the  same  time.  In  all  cases  the  waves 

at  the  response  frequency  and  Its  1/2  subharmonic  become  largely  extinct  upon 
completion  of  the  merging,  while  the  wave  at  the  forcing  frequency  becomes 
distinctly  dominant,  reaching  saturation  at  a  self-limiting  amplitude  consistent 
with  its  wavelength.  This  saturated  wave  combines  with  the  attendant  commen¬ 
surate  distortion  of  the  mean  flow  to  yield  a  flow  structure  which,  when  viewed 
in  coordinates  scaled  to  the  local  mixing  layer  thickness,  is  largely  identi¬ 
cal  I  to  those  encountered  at  subharmonic  saturation  in  modes  1  and  2.  Self-pre¬ 
serving  development  of  the  flew  is  then  regained  downstream  of  the  station 
where  the  vortex  triad  coalesces,  in  agreement  with  the  observations  of  Ho  and 
Huang  (1980). 

The  response  in  mode  4  is  qualitatively  similar  to  that  in  mode  3  to  the  extent 
that  saturation  of  the  response  frequency  is  followed  first  by  amplification  of 
its  1/2  subharmonic  and  then  by  a  three-wave  interaction  with  the  forced  1/4 
subharmonic, which  ultimately  becomes  dominant  and  saturates  at  its  self-limiting 
amplitude.  Self-preserving  development  of  the  flow  then  ensues  downstream  of 
the  station  where  quadruple  vortex  coalescence  is  completed,  again  in  agreement 
with  the  observations  of  Ho  and  Huang  (1980).  The  nonlinear  viewpoint  and 
mechanism  advanced  here  thus  allow  a  unified  interpretation  and  description  of 
the  processes  whereby  a  mixing  layer  subject  to  diverse  initial  conditions  and 
external  forcing  relaxes  to  a  self-preserving  behavior.  The  comments  below  seek 
to  ascertain  the  role  of  that  mechanism  in  the  fully  developed  turbulent  regime. 

In  a  fully  developed  turbulent  flow,  the  dominant,  large  scale  instabilities 
considered  here  are  subject  to  distinct  amplitude  and  frequency  modulations, 
induced  either  by  externally  imposed  disturbances  or  by  secondary  three-dimen-: 
sional  instabilities  intrinsic  to  the  vorticity  patterns  which  evolve  with  the 
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large  scale  structures.  Given  the  nonlinear  nature  of  the  instability  evolution 
process,  these  modulations  result  in  the  seemingly  random  occurrence  In  time 
and  space  of  the  major  observable  events,  e.g.,  vortex  pairings.  Thus,  for 
example,  In  a  naturally  transitioning  mixing  layer,  where  the  initial  amplitude 
and  frequency  of  the  selected  response  disturbance  are  subject  to  modulation, 
the  first  vortex  pairing  occurs  randomly  in  time  and  space  mainly  because  the 
strength  and  separation  of  the  participating  vortices  vary  with  the  modula¬ 
tion  of  the  generating  disturbance.  As  a  result,  the  mean  vorticity  thickness  6^ 
displays  linear  growth  beginning  at/about  the  station  where  the  random 
pairing  occurs  in  the  mean. 

In  a  forced  mixing  layer,  where  the  initial  amplitude  and  frequency  of  the 
response  disturbance  are  determined  uniquely,  the  position  of  first  vortex  coa¬ 
lescence  (either  binary  or  multiple  depending  on  the  response  mode)  becomes 
localized.  Up  to  that  position,  the  thickness  6^  grows  in  the  stepwise 
fashion  predicted,  for  example,  by  the  strictly  twodimensional,  deterministic 
calculations  of  Aref  and  Siggia  (1980)  and  this  paper.  Beyond  that  position 
however,  the  thickness  6^  grows  linearly  because  secondary  instabilities, 
intrinsic  to  the  vorticity  patterns  generated  by  the  first  coalescence,  modu¬ 
late  the  amplification  of  the  subharmonic  responsible  for  the  second  coales¬ 
cence  and  cause  that  event  to  occur  randomly  in  space  and  time. 

Prominent  among  the  modulating  disturbances  are  the  secondary  threedimensional 
instabilities  which  cause  the  formation  of  longitudinal  streaks  (Miksad  1972, 
Brown  and  Roshko  1974),  and  the  spanwise  redistribution  of  the  dominant  mode 
energy  alluded  to  in  Section  5.  Their  effect  is  to  impart  spanwise-variable 
rates/stages  of  evolution/merging  upon  the  vortical  structures  associated 
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with  the  dominant  mode  at  a  ground-fixed  observation  station.  The  Immediate 
consequences  are  twofold:  I)  the  evolving  structures  attain  spanwise 
'  skewed/ branched  configurations  as  observed  by  Browand  and  Troutt  ( 1980) ; 

II)  sequential  structure  realizations  at  a  fixed  observation  station  exhibit 
jitter  in  their  phase  and  amplitude/stage  of  evolution  according  to  some 
probability  density  function  (p.d.f.).  In  the  present  view  the  scale  relation¬ 
ship,  as  well  as  the  mode  and  extent  of  interaction  between  the  dominant  stream- 
wise  and  spanwise  instabilities,  are  deterministic;  randomness  is  relegated 
to  the  spanwise  phase  exhibited  by  the  latter  instability  in  sequential 
realizations  of  the  structures  at  the  observation  station.  Specifically,  the 
spanwise  instability  embodies  the  selective  response  of  the  supporting 
vorticity  concentration  to  the  random  forcing  provided  by  the  prevailing 
three-dimensional  disturbance  environment.  The  receptivity  and  amplification 
characteristics  of  the  vortical  structure  with  respect  to  spanwise  disturbances 
act  as  band  pass  filters  to  select  the  "response".  As  a  result,  the  wave¬ 
length  of  the  spanwise  instability  necessarily  scales  with  that  of  the  locally 
dominant  streamwise  mode  and  with  the  related  local  mean  vortical  thickness  6 

US 

of  the  mixing  layer.  Upon  selection,  the  spanwise  instability  grows  to  an 
amplitude  limited  by  the  nonlinear  requirement  that  the  ensuing  threedimensional 
stretching  of  the  supporting  vortical  structure  remain  bounded.  As  a  result, 
the  p.d.f.  of  the  structures  observed  at  a  fixed  station  possesses  a  well- 
defined  finite  spread;  also,  the  structures  persist  and  cyclically  reproduce 
themselves  with  quasi -two  dimensional  form  at  large  distances  downstream. 

Finally,  if  the  structure  generation  mechanism  and  life  cycle  are  self-pre¬ 
serving  as  proposed  in  this  paper,  so  must  be  their  interactions  with  the 
spanwise  instabilities  and,  therefore,  the  statistics  of  structures  observed 
at  different  uownstream  stations.  The  experimental  findings  that  the  p.d.f. 
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of  vortical  structure  lifetime  In  a  mixing  layer  (Roshko  1976)  as  well  as 
the  spanwlse  distribution  of  cross-correlation  coefficient  (Browand  and 
Troutt  1980)  are  In  fact  self-preserving,  support  this  conclusion.  Thus, 
deterministic  nonlinear  mechanisms  and  bounds  appear  to  control  not  only 
the  twodimensional  evolution  of  the  structures  but  also  their  concurrent 
seemingly  random,  finite  amplitude,  threedimensional  modulations  and  distor¬ 
tions,  which  are  prerequisite  to  the  aforenoted  onset  of  self-preserving 
behavior,  with  linearly  growing  mean  thickness,  following  the  first 
coalescence  of  vortices  In  forced  as  well  as  unforced  mixing  layers. 

A  consequence  of  the  interaction  between  large  scale  streamwlse  and  span- 
wise  instabilities,  and  the  attendant  jitter  of  the  structures,  is  the 
broadening  of  the  low  frequency  peak  In  the  spectra  of  fixed  point  measure¬ 
ments  In  a  turbulent  mixing  layer.  Unfortunately,  this  broadening  does 
not  necessarily  control  the  shape  of  the  low  frequency  spectrum.  The  large 
scale  nonlinear  waves,  which  dominate  the  flow  according  to  the  present  model, 
have  nonsinusoidal  shape  (see  Section  4)  further  complicated  by  the  modu¬ 
lations  attendant  to  the  jitter.  The  coherence  of  such  waveforms  Is  easMy 
lost  in  a  spectral  representation,  which  disregards  phase  information  and 
provides  only  an  averaged  linear  Fourier  decomposition  of  the  signal.  For 
the  signal  shapes  on  hand  the  Fourier  decomposition  typically  includes 
several  artifacts  which  further  broaden  the  spectrum  (Lake  and  Yuen  1978). 
Thus,  the  proposed  dominant  wave  model  of  mixing  layer  dynamics  Is  not 
contradicted  by  the  absence  of  spectral  peaks  in  the  large  body  of  accepted 
measurements  within  turbulent  mixing  layers. 
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An  additional  noteworthy  aspect  of  the  interaction  between  the  large  scale 
streamwlse  and  spanwise  instabilities  in  a  mixing  layer  becomes  manifest 
upon  a  comparison  of  the  mean  thickness  growth  rates  measured  in  experi¬ 
ments  with  those  predicted  by  deterministic,  strictly  two-dimensional 
models,  which  Ignore  the  spanwise  instability  (e.g.  Aref  and  Siggia  1980  \ 

as  well  as  Section  5  of  this  paper).  Remarkably,  theory  and  observation 
are  in  reasonable  agreement,  thus  indicating  a  null| average  effect  of  the 
interaction  upon  the  development  jof  the  flow.  Several  additional  infer¬ 
ences  may  then  be  advanced.  First  of  all,  the  twodimensional  models 
adequately  describe  the  distribution  of  Reynolds  stress  in  the  flow.  Con¬ 
versely,  Reynolds  stress  distributions  so  predicted  provide  a  first  principle 
alternative  and/or  complement  to  flow  analysis  by  the  generally  accepted 
second  order  closure  models,  whose  success  often  depends  on  the  suitable 
empirical  selection  of  third  order  correlation  coefficients.  Secondly,  the 
cascade  of  energy  in  threedimensional  wavenumber  space  provides  an  important 
contribution  to  the  level  and  distribution  of  turbulent  kinetic  energy  in 
the  flow,  but  has  little  Influence,  either  qualitative  or  quantitative,  on 
the  evolution  of  the  large  scale  structures  and  on  the  distribution  of 
Reynolds  stress.  Since  the  distributions  of  all  second  order  turbulence 
quantities  in  a  mixing  layer  become  self-preserving  after  the  second  vortex 

I 

pairing  (Browand  and  Latigo  1379)  »  the  cascade  must  be  intrinsic  to,  and 
bound  with,  the  evolving  vorticity  patterns  of  the  large  scale  structures, 
and  must  evolve  subject  to  the  same  broad  constraints  previously  indicated 

*  We  conjecture  that  the  low  frequency  and  inertial  subranges  of  spectra 
become  also  self  preserving  after  the  second  vortex  pairing. 
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for  the  large  scale  spanwise  Instability  —  i.e.,  wave  vectors  selected 
by  the  receptivity  and  amplification  characteristics  of  the  supporting 
advected  flow  patterns,  phase  determined  case  by  case  by  the  prevailing 
random  disturbance  environment,  amplification  limited  by  nonlinear  bounds 
on  the  ensuing  distortion  of  the  supporting  flow  patterns,  effect  confined 
to  localized  distortions  of  the  dominant  carrier  waveshape  and  the  asso¬ 
ciated  distributions  of  fluctuations.  The  overall  energy  of  the  carrier 
wave  and  its  evolutionary  interaction  with  the  mean  flow  are  then  influ¬ 
enced  only  slightly,  and  the  self-preserving  statistics  of  the  cascade 
process,  as  well  as  its  negligible  effect  upon  the  dynamics  of  the  large 
scale  structures  and  the  development  of  the  flow,  are  rationalized. 

The  arguments  above  support  the  premise  of  the  proposed  model  that  the. large 
scale  structures  and  the  mean  flow  of  a  mixing  layer  constitute  a  largely 
i  self-contained  nonlinear  twodimensional  flow  system,  which  evolves  by  a 

self-reproducing  deterministic  limit  cycle.  The  cycle,  although  subject 
to  modulation  by  intrinsic  self-limiting  instabilities,  becomes  rapidly  ob¬ 
livious  of  initial  conditions.  Control  of  the  flow  then  seems  possible 
only  by  repeated  local  interference  with  either  the  mean  flow  or  the  domi¬ 
nant  modes  associated  with  the  large  scale  structures.  For  example,  one 
anticipates  that,  in  line  with  the  findings  of  Ho  and  Huang  (1980),  the 
,  growth  rate  of  the  mean  mixing  layer  thickness  at  a  selected  downstream  dis¬ 

tance  is  enhanced  by  in  phase,  small  amplitude,  external,  forcing  of  either 
the  locally  dominant  wave  frequency  or  its  nondispersive  fractional  subharmonics 
Conversely,  one  expects  that  the  growth  rate  is  diminished  if  the  frequency 

1 

of  the  forcing  is  preserved,  but  its  phase  is  reversed.  Receptivity  consid¬ 
erations  together  with  the  proposed  model,  and/or  slight  extensions  thereof, 

I 

1 


NYU/DAS  81-06 


should  provide  a  realistic  tool  for  evaluating  the  effectiveness  of  such 
control  concepts  and  for  designing  attendant  demonstration  experiments. 
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FIGURE  1:  The  preexponential  factor  A.(t)  and  the  amplification  rate 

3-(t)  of  the  fundamental  wave  according  to  the  extended  Meksyn 
aAd  Stuart  model. 
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FIGURE  2: 
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FIGURE  3: 


Experimentally  observed  downstream  growth  of  the  Integrated  u' 
energy  of  disturbances  in  a  mixing  layer  forced  at  the  dlmen-  ’ 
sionless  frequency  8  =  0.222.  x,  0  ;  o,  20;  a,  38;  P,( 1/2)8 
7, (3/2)0.  From  Mi ksad  (1972). 
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FIGURE  9-.  The  evolving  normalized  stream  function  C,(y,t)  of  the  fundamental 
wave  motion  at  selected  times  following  its  saturation.  Nonlinear 

model  with  subharmonics  suppressed,  oo  o  t  =  13.5*,  - - - , 

t  *  19.15; - ,  t  =  22.15;  —  - - -  t  *  25.75;  — -  t  =  30.55. 
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FIGURE  10:  The  evolution  of  the  mean  velocity  profile  U(y,t)  at  selected 
times  following  excitation  of  the  1/2  subharmonic.  Nonlinear 
model,  o  o  o,  t.=  17;  -. - ,  t.  =  .22,4;  — - .  t.  =  30.2. 


FIGURE  12:  Time  evolution  of  the  normalized  integrated  energy  E  (t)  of  the 
fundamental  (n  =  1)  and  of  the  subharmonic  (n  =  l/2)^nthe  latter 
being  excited  at  t  =  17.  Nonlinear  model. 


